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Abstract 


We  present  an  evaluation  of  the  potential  improvements  in  seismic  disturbance  rejection  to 
be  obtained  by  using  active  variable  damping  control  in  a  structure.  Using  the  responses  to 
seismic  excitation  of  an  optimally  controlled  variable  structure  and  of  a  minimax  optimally 
designed  fixed  structure,  we  obtain  an  upper  bound  on  the  achievable  performance  and  a  lower 
botmd  on  the  acceptability  of  a  control  system  for  a  variable  damping  structure. 

Our  numerical  experiments  suggest  the  following  three  conclusions: 

(i)  A  minimax  optimal  designed  stmcture  gives  very  good  seismic  disturbance  suppression, 
not  only  for  the  earthquakes  used  in  its  design,  but  also,  for  other  earthquakes  of  similar  intensity. 

(ii)  The  use  of  variable  structure  control  is  not  likely  to  result  in  smaller  maximum  interstory 
drifts  than  exhibited  by  a  fixed  minimax  designed  structure.  However,  the  use  of  variable  struc¬ 
ture  control  may  be  somewhat  more  effective  in  in  the  reduction  of  peak  and  average  accelera¬ 
tions.  Overall,  controlled  variable  structures  are  likely  to  perform  best  when  the  earthquakes  are 
moderate  to  severe,  and  at  sites,  such  as  landfills  and  dry  lake  beds,  where  resonances  can  be 
expected,  but  the  resonance  frequency  cannot  be  estimated  in  advance. 

(Hi)  The  gap  between  the  upper  and  lower  bounds  is  rather  small,  which  makes  designing  a 
feedback  law  that  results  in  performance  superior  to  that  of  a  minimax  designed  structure,  very 
difficult.  The  best  choice  appears  to  be  in  the  form  of  a  continuous  moving  horizon  control  law, 
using  a  horizon  of  under  0.2  seconds.  An  implementation  of  such  a  control  law  will  require 
determination  whether  local  ground  motion  can  be  predicted  up  to  0.2  sec.  ahead,  using  ground 
motion  monitoring  sensors  located  a  small  distance  away  from  the  site. 
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1.  introduction 


The  problem  of  controlling  seismically  excited  vibrations  in  a  terrestrial  stmcture 
has  features  that  make  it  significantly  different  from  the  problem  of  controlling  vibrations 
in  a  structure  caused  by  wind  loading,  or  traffic,  or  from  of  controlling  vibrations  in  a 
flexible  space  structure  caused  by  docking  maneuvers  or  space  debris  impacts.  Seismic 
disturbances  are  of  short  duration,  they  are  potentially  of  destructive  intensity,  and  their 
occurrences  are  separated  by  long  periods  of  quiescence.  In  addition,  since  seismic  load¬ 
ing  is  applied  only  to  the  base  of  a  structure,  many  of  the  methods  of  generating  control 
forces  for  the  elimination  of  wind  loading  and  other  t)q)es  of  excitation  are  not  very 
effective  in  dealing  with  seismic  disturbances. 

The  use  of  passive/active  control  to  reduce  the  damage  caused  by  earthquakes  to 
buildings  and  other  structures  has  become  an  area  of  considerable  theoretical  interest. 
There  are  several  recent  collections  of  papers  on  structural  control  [Lei.l,  Lei.2,  Cho.l, 
Wen.l]  and  many  articles  containing  descriptions  of  control  techniques  applicable  to  the 
seismic  damage  control  problem.  For  a  representative  sample  of  papers  dealing  with  the 
control  of  structures  subject  to  wind,  traffic,  and  earthquake  disturbances,  see,  e.g., 
[Abd.l,  Cha.l,  Kel.l,  Mar.l,  Kob.1,2,  Mas.l,  Roo.l,  Yan.2  Yao.l],  as  well  as  the  over¬ 
view  paper  by  T.  T.  Soong  et  al  [Soo.2].  The  recent  text  by  T.  T.  Soong  [Soo.l]  demon¬ 
strates  that  the  field  of  structural  control  is  reaching  a  certain  level  of  maturity.  An 
interesting  aspect  of  the  research  on  the  control  of  terrestrial  structures  is  the  range  of 
actuation  options  that  are  being  considered,  such  as  active  base  isolation  systems  (see, 
e.g.,  [Kel.l]),  the  use  of  active  tendons  (see,  e.g.  [Roo.l]),  active  mass  damping  systems 
(see  e.g.,  [Cha.l,  Kob.2]),  semi-active  impact  damping  systems  (see,  e.g.  [Deh.l])  and 
active  variable  structure  control  [Kob.l]. 

The  control  laws  used  in  structural  control  often  make  use  of  state  feedback  in  the 
form  of  instantaneous  optimal  control  laws,  based  on  the  minimization  of  linear  quadratic 
cost  functions  (see,  e.g.  [Yan.3,  Soo.l,  Che.l,  Sat.l]),  independent  modal  control,  (see, 
e.g.  [Mei.l,  Yan.4,  Mar.l]),  pole  assignment  (see,  e.g.  [Abd.2]),  sliding  regimes  (see,  e.g. 
[Wan.l]),  bang-bang  control  (see,  e.g.  [Kaw.l]),  Lyapunov  function  based  methods  (see, 
e.g.  [Lee.l,  Kel.l]),  and  pulse  control  (see,  e.g.  [Deh.l,  Pru.l,  Mas.l]). 

For  reasons  of  reliability,  the  energy  for  operating  the  control  system  of  a  seismic 
resistant  structure,  ought  to  be  stored  locally.  The  actuator  choices  have  a  considerable 
impact  both  on  the  energy  consumption  and  on  the  potential  effectiveness  of  the  resulting 
control  system.  For  many  active  control  schemes,  such  as  active  base  isolation  and 
active  tendons,  because  of  the  enormous  weight  of  a  multistoried  structure,  the  energy 
requirements  can  be  very  large  indeed,  and  hence  the  development  of  safe  energy  storage 
methods  for  these  types  of  actuators  remains  an  open  problem.  The  actuators  proposed  in 
[Deh.l,  Kob.l]  are  distinguished  by  the  fact  that  they  use  very  little  energy,  which  is 
highly  desirable,  but  they  raise  the  question  as  to  how  effectively  one  can  control  a  struc¬ 
ture  without  applying  considerable  external  forces. 
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Since  in  a  variable  structure  one  cannot  apply  external  control  forces,  the  function 
of  a  variable  structure  control  system  is  to  produce  ah  impedance  mismatch  between  the 
ground  motion  and  the  structure,  so  as  to  minimize  the  amount  of  ground  motion  energy 
that  gets  pumped  into  the  stincture.  The  success  of  such  a  control  law,  therefore,  is 
bound  to  depend  considerably  on  the  spectral  properties  of  the  ground  motion.  The  pur¬ 
pose  of  this  study  is  to  establish  the  capabilities  of  conttol  laws  for  variable  damping 
structures  as  well  as  for  the  simple  variable  stiffness  structure  proposed  in  [Kob.l]. 

In  section  2.1,  we  present  a  model  of  the  planar  variable  damping  stracture  that  we 
wish  to  study.  Since  the  differential  equations  describing  this  model  are  stiff,  simulations 
using  this  model  are  rather  time  consuming.  We  therefore  proceed  to  develop  a 
simplified  model  that  may  be  useful  for  use  in  various  preliminary  studies,  involving  the 
use  of  sampled-data  control  laws.  In  Section  2.2,  we  examine  the  behavior  of  the  most 
drastic  simplification  of  our  model,  in  which  the  mass  of  the  auxiUaiy  structure  is 
neglected,  and  the  variable  dampers  are  assumed  to  function  as  on-off  clamps  only.  In 
the  process  we  examine  the  behavior  of  the  original  variable  stiffness  configuration  pro¬ 
posed  in  [Kob.l],  in  which  an  auxiliary  structure  is  connected  to  a  main  stracture  by 
means  of  controllable'  clamps,  i:e,  the  stiffness  variation  is  of  the  on-off  variety.  We 
show  by  analysis  that,  when  operated  in  sampled-data  relay  mode,  when  the  sampling 
period  goes  to  zero,  such  a  structure  can  be  made  to  approximate  a  structure  with  stiff¬ 
ness  values  that  lie  anywhere  in  the  convex  hull  of  those  of  the  main  stracture  and  the 
all-clamped  stracture.  However,  our  experimental  studies  show  that  a  variable  stracture, 
that  is  operated  in  sampled-data  relay  mode,  exhibits  very  sharp,  short  duration  accelera¬ 
tion  peaks.  Since  these  may  cause  buckling,  the  use  of  such  a  control  method  in  practice 
may  be  problematic.  We  therefore  abandoned  this  model,  and,  in  Section  2.3,  we  con¬ 
struct  a  somewhat  more  sophisticated  simplified  model,  in  which  the  mass  of  the  auxili¬ 
ary  structure  is  still  neglected,  so  as  to  avoid  creating  stiff  differential  equations,  but  the 
on-off  clamp  action  of  the  variable  dampers  is  replaced  by  a  smoothed  out  action.  This 
model  turned  out  to  be  adequate  for  use  in  preliminary  studies,  but  not  good  enough  to 
replace  the  original  model  entirely. 

In  Section  3,  we  establish  performance  bounds  for  a  controlled  variable  damping 
structure.  We  begin  in  Section  3.1  by  chosing  an  energy-based  performance  criterion 
that  is  close  in  spirit  to  the  quadratic  criteria  used  in  many  feedback  laws.  In  Section  3.2, 
we  then  establish  upper  bounds  on  achievable  performance  for  five  well  known  earth¬ 
quakes  by  determining  the  response  of  a  variable  damping  structure  with  optimal  damp¬ 
ing  modulation  (using  an  algorithm  from  [Bak.l]),  that  is  computed  under  the  unrealistic 
assumption  that  the  ground  motion  is  known  in  advance.  Then,  in  Section  3.3,  we  use  an 
optimization  algorithm  (see  [Pol.  1,2])  to  compute  fixed  damping  coefficients  that  are 
minimax  optimal,  for  our  energy  criterion,  with  respect  to  four  of  the  earthquakes  con¬ 
sidered.  The  fifth,  the  El  Centro  1940,  earthquake  was  left  out  to  be  used  as  a  control 
case,  to  determine  whether  a  stracture  that  has  been  mininjax  tuned  for  a  given  set  of 
earthquakes  performs  well  on  other  earthquakes.  The  performance  of  this  minimax 
designed  fixed  stracture  determines  a  lower  acceptable  performance  bound. 
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Since  the  performance  gap  between  the  unrealistic  optimally  controlled  structure 
and  the  minimax  designed  stmcture  was  found  to  be  small,  we  decided  to  explore  only 
the  behavior  of  moving  horizon  control  laws  (see,  e.g.  [Pol.3]),  since  they  have  the  best 
chance  of  giving  behavior  that  falls  in  the  gap  between  our  upper  and  lower  performance 
bounds.  In  Section  4,  we  present  our  moving  horizon  control  laws:  the  first  is  of  the 
sampled-data  type,  while  the  second  one  is  a  continuous  control  law. 

In  Section  5  we  evaluate  our  numerical  results.  First,  we  found  that  the  correlation 
between  our  energy  criterion  and  such  basic  requirements  as  keeping  interstory  drifts 
within  prescribed  bounds  and  floor  accelerations  small,  was  adequate,  though  by  no 
means  exact  Second,  we  found  that  the  minimax  designed  structure  performed  very 
well,  even  on  the  El  Centro  1940  earthquake  that  was  not  used  in  its  computations. 
Finally,  we  found  that  the  use  of  a  continuous  moving  horizon  control  law,  with  an  accu¬ 
rate,  0.2  sec.  prediction  of  the  ground  motion  is  capable  of  yielding  performance  that  is 
better  than  our  minimax  design. 

Our  concluding  thoughts  are  presented  in  Section  6.  In  particular,  we  arrive  at  two 
conclusions.  The  first  is  that  sophisticated  minimax  design  (see,  e.g.,  [Pol.  1,2]),  taking  a 
fairly  small  number  of  ground  motions  into  account,  is  likely  to  yield  structures  that 
withstand  seismic  disturbances  much  better  than  structures  designed  by  standard  tech¬ 
niques.  The  second  conclusion  is  that  controlled  variable  structures  are  likely  to  be  per¬ 
form  best  at  sites,  such  as  landfills  and  dry  lake  beds,  where  resonances  can  be  expected, 
but  the  resonance  frequency  cannot  be  estimated  in  advance. 

2.  System  Models  for  Numerical  Experimentation 


As  we  will  see,  the  natural  model  for  a  variable  structure  is  numerically  stiff,  and 
hence  numerical  experiments  using  this  model  are  rather  time  consuming.  We  will  there¬ 
fore  develop  a  simplified  model  which  results  in  much  reduced  simulation  times  and 
which  is  adequate  for  the  tedious  preliminary  numerical  experiments  that  we  need  to  per¬ 
form  in  order  to  establish  effective  weights  in  our  cost  function.  As  we  will  see,  the  sim¬ 
plest  model  that  one  can  propose  has  nice  theoretical  properties,  but  does  not  sufficiently 
accurately  model  the  acceleration  response.  Hence,  a  somewhat  more  complex,  but  still 
numerically  well  conditioned,  simplified  model  had  to  to  be  adopted. 

2. 1.  Variable  Structure  System  Model 

We  will  consider  the  control  of  a  planar  variable  three  story  structure  consisting  of  a 
main  structure  and  an  auxiliary  structure  that  are  linked  through  variable  damping  ele¬ 
ments,  as  shown  in  Fig.  2.1.1.  We  will  assume  that  the  mass  of  the  main  structure  is  con¬ 
centrated  in  its  floors,  and  that  the  mass  of  the  auxiliary  structure  is  concentrated  at  the 
points  of  connection  of  the  auxiliary  structure  to  the  variable  dampers.  Hence  our 
mathematical  model  corresponds  to  the  situation  in  Fig.  2.1.2. 
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Fig.  2.1.1.  Configuration  of  Variable  Damping-Stiffness  Structure 


Fig.  2.1.2.  Idealization  of  Variable  Damping-Stiffness  Structure 


We  will  model  the  variable  structure,  as  a  planar  structure  with  horizontal  displace¬ 
ments  only,  in  terms  of  absolute  coordinates,  as  follows: 

Mm  +  Ci(f)^Kx(t)+  i  =  (2.1.1a) 

i=l 

where  the  vector  a:  e  R®  is  given  by 

x{t)  =  {xXt),x^{t),x^{t),x}(t),x^{t),x^(t)),  (2.1.1b) 

whCTe  the  xj(r),  i  =  1,2,3,  are  the  absolute  displacement  of  the  floors  of  the  main 
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Structure,  the  (t),i  =  1 , 2, 3,  are  the  absolute  displacement  of  the  masses  of  the  auxili¬ 
ary  structure,  M  is  the  mass  matrix,  C  is  the  fixed  damping  matrix,  K  is  the  stiffness 
matrix,  p'  (t),i  =  1, 2, 3,  are  the  variable  damper  coefficients,  11^,1  =  1, 2, 3,  are  location 
matrices  for  the  variable  dampers,  Jg  is  a  6  x  1  vector  with  each  element  equal  to  1,  and 
Xp(0  and  are  ground  motion  and  ground  velocity  obtained  from  an  earthquake 

o  o 

record.  We  assume  that  at  rest  the  positions  of  the  floors  and  of  the  auxiliary  masses  are 
zero,  in  absolute  coordinates. 

The  parameter  values  for  our  model  are  given  in  Table  2.1.1,  where  /  =  1, 2, 3. 

Next,  defining  the  state  vector  z  e  by  z(t)  =  (x(t),x(t)),  and  e(t)^  R^  by 
e  (t )  =  (Xg  it),Xg  (t)),  equation  (2.1.1)  can  be  rewritten  in  the  canonical  state  space  form 

3  . 

z{t)  =  Az(t)+'^  u^{t)BiZ{t)+De{t).  (2.1.2) 

i=l 

This  system  of  differential  equations  turns  out  to  be  rather  stiff,  particularly  when 
the  variable  damping  coefficients  are  high,  because  the  mass  of  the  auxiliary  structure  is 
much  smaller  than  the  mass  of  the  main  structure.  Hence  numerical  simulation  times 
tend  to  be  rather  high.  Therefore,  in  the  next  two  subsections,  we  will  attempt  to  develop 
a  simplified,  numerically  well  conditioned  model  that  can  be  used  in  early  stages  of  con¬ 
trol  system  design. 


2.0  X  10^  kg 

K 

1.2  X  10^  kg 

1.2  X  10^  kg 

K 

1.96X  lO^V/m 

1.96xlO''V/m 

2.45X  lO^N/m 

V^s 

0.004/i:i  Ns/m 

Mi 

l.OxlO^^g 

K 

9.8xl0''iV/m 

0.004Kj,Ns/m 

Table  2.1.1.  Parameter  Values  for  Variable  Stmcture. 
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2.2.  Simplified  Variable  Structure  System  Model  I 

Since  one  of  the  control  laws  that  we  propose  to  explore  is  a  moving  horizon 
sampled-data  control  law,  we  begin  by  exploring  a  simplification  of  the  model  (2.1.1),  in 
which  the  mass  of  the  auxiliary  structure  is  neglected,  and  the  variable  dampers  are 
modeled  as  on-off  clamps.  We  will  present  a  theoretical  analysis  that  displays  the  capa¬ 
bilities  of  the  resulting  variable  structure  with  a  sampled-data  clamp  switching  law,  as 
well  as  a  set  of  numerical  experiments,  designed  to  evaluate  the  validity  of  this  simplified 
model  of  the  dynamics  and  control  action. 

Consider  the  simple  structure  in  Fig.  2.1.1,  and  suppose  that  the  variable  dampers 
act  as  on-off  clamps.  Since  each  clamp  can  be  either  open  or  closed,  at  each  point  of 
time,  the  behavior  of  the  structure  is  determined  by  one  of  the  eight  possible  combina¬ 
tions  of  the  clamp  states.  To  keep  track  of  this  situation,  we  introduce  the  clamping  vec¬ 
tor  c(r)  e  r  ^  0,  with  die  components  c'(/),  i  =  1,2, 3,  assuming  only  the  values  of 
0  or  1,  i.e.,  c  (/)  e  C,  where 


'o' 

fol 

fo] 

fol 

fil 

fil 

fil 

fil 

■> 

0 

1 

0 

1 

0 

1 

0 

1 

•*  . 

(2.2.1) 

loj 

LqJ 

UJ 

LiJ 

Iq) 

0 

V  ^ 

UJ 

LiJ 

When  c‘(t)=  1,  the  i-th  clamp  is  closed,  otherwise  it  is  open.  We  will  assume  that  c(t) 
can  only  change  at  the  sampling  times  iT,  k  =  0, 1,2,3,. . .,  where  T  is  the  sampling 
period. 

Hence,  using  coordinates  relative  to  the  inertial  frame  (i.e.,  absolute  coordinates), 
for  r  >  0,  to  indicate  the  dependence  of  the  motion  of  the  system  on  the  clamping  vector 
c(r),  we  will  denote  the  resulting  position  of  the  i-th  mass  of  the  main  structure  by 
xj(r,c),i  =  1,2,3,  and  the  position  of  the  f-th  clamp  by  Xa(f,c),  I  =  1,2,3.  Finally,let 
x/(r,c)  =  x^(r),  the  absolute  position  of  the  ground.  We  assume  that  the  absolute 
groimd  motion  is  caused  by  a  seismic  disturbance,  and  that  it  is  at  least  twice  continu¬ 
ously  differentiable. 

Next,  for  i  =  1, 2, 3,  let 


Axj(r ,  c )  ^  xj(r ,  c )  -  ,  c ), 

(2.2.2a) 

Axj(t ,  c )  4  xj(r ,  c )  -  c ), 

(2.2.2b) 

(2.2.2c) 

c\t)  =  Ax>,c)  =  Ax,%.c)  =  A;c>,c)  =  0, 

(2.2.2d) 

(2.2.2e) 

where  the  last  two  relations  are  introduced  so  as  to  make  it  possible  for  us  to  write  down 
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the  equations  of  motion  of  the  controlled  structure  in  compact  form,  as  follows, 
main  structure.  For  /  =  1, 2, 3, 

(r ,  c ) + (nj  +  c  *  (t  )|xi  )Axj  (r ,  c ) + /sTj  Ax/(r ,  c ) + c  ‘  (r  )/s:^  Ax^  (r ,  c ) 


+K^-^Axi-\t,c )  +  c^-\t  )Kir^Axi-\t,c ), 

(2.2.3a) 

auxiliary  structure.  For  i  =  1, 2, 3, 

Xait,c)  =  C^it)x^it,C)  +  il-c'it))[Xs  it,c)  Ax^(f,c)]. 

(2.2.3b) 

We  assume  that 

xi(0)  =  ii(0)  =  0,  i  =  1,2, 3, 4, 

(2.2.3c) 

xi(0)=xi(0)  =  0,/=  1,2,3, 

(2.2.3d) 

and  that  the  ground  motion  (Xg  (t),  Xg  (t))  is  given. 

If  we  define  the  state  vector  z  e  by 

z(0  4  (2.2.3e) 

and  define  the  ground  excitation  function  e(t)  =  (Xg(t),Xg(t)),  then  we  conclude  from 
(2.2.2a),  (2.2.2b)  that  z(r)  satisfies  a  differential  equation  of  the  form 

z(t)  =  A  (c  (t))z  it)  +  Bic  it))e  it),  (2.2.4a) 

where  the  matrices  A  (c  (r ))  and  B  (c  it ))  have  the  form 

3  . 

A  (c  (r ))  =  A  0  +  X  c  ‘  (OA,- ,  (2.2.4b) 

i=l 

Bicit))  =  BQ  +  c\t)B^,  (2.2.4c) 

where  the  matrices  A,-,  i  =  0, 1,2, 3,  and  Bq,  5 3  are  constant.  Note  that  although  the 
clamp  status  c(r)  is  not  necessarily  continuous,  by  the  fundamental  theorem  of  differen¬ 
tial  equations,  the  state  trajectory  z  (t )  is  continuous. 

To  establish  the  potential  capability  of  a  variable  damping-stiffness  strucmre,  we 
will  show  that  the  use  of  a  relay  clamp  switching  law  enables  us  to  endow  the  structure 
with  any  dynamics  in  a  certain  set.  For  this  purpose,  we  assume  that  we  need  to  analyze 
the  behavior  of  the  structure  over  the  time  interval  [0,7^],  and  that  our  relay  clamp 
switching  law  is  based  on  a  sampling  time  T  4  T^INM,  where  N  and  M  are  positive 
integers  such  that  [0,7^]  can  by  divided  into  N  intervals  4  =  [kMT ,ik+l)MT),  for 
k  =  0, 1,...,N-1.  The  clamping  function  c(t)  then  has  the  property  that  each  of  the 
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components  c‘(t)  is  constant  over  each  sample  period  [kT,ik+l)T)  where 
^  e  { 0, 1, . .  .,NM-1 } .  We  then  define  the  modulation  function  m ;  R  ->  [0, 1]^,  by 

m^(t)=pilM,teIk,  i  =  1,2,3,  ^  =  0, 1,.. (2.2.5) 

where  0<p^<M  is  the  number  of  subintervals  in  4  for  which  the  clamp  connected  to 
floor  i  is  closed,  (ie.,  c‘(t)=  1). 

Now  consider  the  “averaged”  dynamics: 

(0  =  ^av  (t)z^(t)  (t)e  (t),  (2.2.6a) 

where 

3 

^av (^)  “'^0'^  £  ni  (t )A^ ,  (2.2.6b) 

1=1 

Bav(t)  =  BQ  +  m\t)B^  .  (2.2.6c) 

Let  6A,(0^[c‘(t)-m‘(t)M/,  /  =  1,2,3,  and  6fi3(t) ^ [c^(t)-m^(0]fi3.  Then, 
assuming  that  z  (0)  =  z^  (0)  =  0,  we  find  that  for  f  e  [0,  Ty  ], 

z (0  - Zflv (0  =  [A  (c (s ))z (s )  -  A^^ (s)z^^ (s)  +  (B  (c (s ))  - B^ (s ))e (s )]ds 

=  l^^av(^)[z(s)~Zav(s)]ds 

+  j^[A(c(s))-A^(s)]z(s)ds 

+  Jo  [B  (c  (s ))  -  (s  y\e  (s)ds  .  (2.2  Ja) 

Now, 

3 

A  (c  (s ))  -  /lav  (^  )  =  Z  S/1,  (s ),  (2.2.7b) 

i=l 

and 

B(c(s))-B^^(s)  =  8B2(s).  (2.2.7c) 

Note  that  for  any  interval  4,  k  =  0, 1, . .  .,N-1,  and  i  =  1,2,3, 

W  is)-m\s)']AiZ(kMr)ds  =  \pIT -Mr(pllM)]AiZ(kMt)  =  0.  (2.2.7d) 

Now,  let  Lj5  be  the  Lipschitz  constant  for  z(  )  on  [0, 7y  ],  and  ilTa  <  <»  be  a  constant  such 
that  1/1, •  I  <  Ka  for  i  =  1, 2, 3.  Then  for  ^  =  0, 1, . . ., N-\,  and  /  =  1, 2, 3, 

*IiLt  ^  (s  life  I  ^  “  z  (kMT)\ds 

+  iJ^^^5A,(^)z(/irMr)d:yi 

-kMTlds  (2.2.1c) 

For  any  t  e  [0,7^^],  there  exists  r  e  N,  0  ^  r  <iV-l  such  that  rMT  ^  t  ^  (r+l)MT  and 
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hence, 

t  '■■"1  3  ,(k+DMT  3  I 

If  [/l(ca))-A„(i))z(s)*<S  2  H(s)2(^)*II+  SA,a)z(s)*l 

k=Oi=\  i=l 

<  3rK^L,M'^'^  +  3(t  -  rMIW^K, 

<  3NK^ 4  M^T'^  +  3MTK^ K, 


=  3MTK^[L,Tf+K,l  (2.2.7f) 

where  <  oo  is  such  that  llz  (r)l  <  for  all  t  e  [0,Tf].  Similarly,  with  4  the  Lipschitz 
constant  for  e(  )  on  [0,Tf],  Kf,  <°°  such  that  and  Kg  <oo  such  that 

lie  (t)i<Kg  for  all  t  e  [0,4 ],  we  get  that 

[B  (c  (s ))  - is )]e  is)dsl< MTKt,  [Lg 7>  +  ATJ.  (2.2.7g) 

Next,  since  is  bounded  for  all  t,  we  can  assume  that  there  exists  a  constant 

Kav<°°  such  that  (5  )l  ^  AT^v  for  all  5  e  [0,7^:],  and  setting 

K  =  3K^  [4  4  +  ATJ  +  ATf,  [4  4  +  ATJ ,  we  conclude  from  (2.2.7a),  (2.2.7f),  and  (2.2.7g) 
that 

hit)-z^,it)i<K^fhis)-z^,is)lds  +KMT,  (2.2.7h) 

It  now  follows  from  the  Bellman-Gronwall  Lemma  that 

T 

hit)-  it)\  <  KMTe^'^'^f  =  A:-^e^‘■''^^  (2.2.7i) 

where  we  have  used  the  fact  that  MNT  =  7y . 

Thus  we  see  that  as  the  number  of  samples  N ,  in  [0,4]*  S^es  to  infinity,  the 
behavior  of  the  relay  clamp  switched  structure  becomes  the  same  as  that  of  the  averaged 
dynamics.  Note  that  the  error  between  the  true  and  averaged  dynamics  does  not  neces¬ 
sarily  go  to  zero  as  the  sample  size,  T,  goes  to  zero.  In  the  case  where  N  remains 
bounded,  we  have  only  shown  a  fixed  upper  bound  on  the  error.  However,  when  both  N 
and  M  go  to  infinity,  in  the  limit,  the  error  goes  to  zero  and  the  modulation  function  m  (t) 
can  attain  any  value  in  a  dense  subset  of  [0, 1]^.  In  the  special  case  where  m  (r)  =  m*  is 
constant,  the  relay  clamp  switched  structure  behaves  as  a  linear  system  with  dynamics 
defined  by  m* . 

numerical  evaluation  of  the  model.  We  carried  out  a  numerical  experiment 
using  the  model  in  this  section  so  as.  to  determine  its  validity.  We  used  a  sampling  time 
T  =  0.2  seconds,  a  final  time  4  “  ^  seconds,  the  earthquake  ground  motion  supplied  by 
the  El  Centro  1940  record,  and  a  command  signal  constracted  using  a  simple  moving 


-9- 


2.2.  Simplified  Variable  Structure  System  Model  I 


horizon  control  law.  Plots  of  the  resulting  accelerations  of  the  floors  of  the  main  struc¬ 
ture,  and  the  clamping  vector  c(t),  are  shown  in  Fig.  2.2.1  and  Fig.  2.2.2.  Note  that  all 
plots  included  in  this  report  follow  the  convention  that  the  solid  line  refers  to  floor  1  (top) 
of  the  structure,  or  to  the  first  component  of  the  vector  being  plotted;  the  dashed  line 
refers  to  floor  2  (middle)  of  the  structure,  or  to  the  second  component  of  the  vector  being 
plotted;  and  the  dotted  line  refers  to  floor  3  (ground)  of  the  structure,  or  to  the  third  com¬ 
ponent  of  the  vector  being  plotted.  Referring  to  Fig.  2.2.1  and  Fig.  2.2.2,  we  see  large, 
short  duration  acceleration  peaks  that  occur  at  the  clamp  opening  and  closing  times.  The 
reason  for  the  acceleration  peaks  is  that  the  opening  and  closing  of  a  clamp  causes  impul¬ 
sive  forces  to  act  on  the  main  and  auxiliary  structures.  Since  large  acceleration  peaks 
may  cause  buckling  failure  in  the  stracture,  we  will  adopt  a  law  for  switching  the 
dampers  that  is  gradual,  for  use  in  sampled-data  relay  control.  This  results  in  a  modified 
structural  model  that  we  will  present  in  the  next  subsection. 

2.3  Simplified  Variabie  Structure  System  Modei  ii 


The  numerical  results  presented  in  the  last  section  indicate  that  even  for  preliminary 
studies,  a  more  sophisticated  model  than  the  one  presented  in  Section  2.2,  for  the  struc¬ 
tural  system,  is  required.  We  will  therefore  replace  the  clamps  by  variable  dampers 
whose  damping  coefficients  can  be  varied  from  O.ON  secim  to 

o 

lijjjax  =  1.96  X  10  N  sec  Im .  As  in  Section  2.2,  we  will  assume  that  the  command  signal 
c  (r )  e  IR^  takes  values  in  the  set  C,  defined  in  (2.2.1),  and  that  it  can  only  change  at  the 
sampling  times  kT,k  =  0, 1, 2, 3, ... ,  where  T  is  the  sampling  period.  Since  the  damp¬ 
ing  coefficients  depend  on  the  command  signal,  we  will  denote  them  by  3, 

and  we  will  assume  that  they  respond  to  the  command  signal,  c  (t ),  as  follows: 


I  MmaxCa -mrrj  fort  e  [kT,kT  +T,] 

IVax  fOTt^(kT+T„ik+l)T), 

if  c‘(kT  -  e)  =  0  and  c*(kT)  =  1  (with  e  e  (0,  T)); 


(2.3.1a) 


P™„(1  -  (r  - kT)/T^Y  for /  e  [kT. kT +T^] 

0  fort  ^  (kT +Tc,(k  +  1)T), 


(2.3.1b) 


ifc‘(/tr-e)=  landc‘(itr)  =  0; 


11^(0  =  0,  for  t  e  [kT,  (k  +  1)7),  (2.3.1c) 

if  c^(kT  -  e)  =  0  and  c‘(kT)  =  0;  and,  finally, 

=  Mmax.  for  t  e  [kT,  (k  +  1)7),  (2.3.1d) 

if  c‘  (kT  -  e)  =  1  and  (kT)  =  1.  The  exponent  r  e  IN  will  be  chosen  on  the  basis  on 
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numerical  experimentation. 

The  introduction  of  a  variable  damper  modifies  the  equations  of  motion  (2.2.3a, b)  as 
follows. 


main  structure.  For  i  =  1, 2, 3, 

it,c)  +  K^Axlit,c)  +  + 


=  K^-^AxtHt  ,c)  +  ,c)  +  ,c)  +  ,  c  ),(2.3.2a) 


i-l  AO‘-l/ 


auxiliary  structure.  For  i  =  1, 2, 3, 

V^iit)(.Xsit,c)-Xa(t,c))  =  (2.3.2b) 

Let  |Xc(0  =  (M-c(0.  vector  of  damping  coefficients  correspond¬ 

ing  to  the  three  clamps.  As  before,  by  defining  the  state  vector  z  e  by 

z(0  ^  (x,\t),x,\t),x^\t),x,\tlx,\t),x^(t),x^{t),x^it),x^it)),  (2.3.2c) 

and  the  ground  motion  function  e(t)  =  (Xg(t),Xg(t)),  then  from  (2.3.2a,b),  we  can  show 
that  z  (r )  satisfies  a  differential  equation  of  the  form 

z(r)  =  A(p(i,c))z(0+^(|Xc(0)e(0,  (2.3.3a) 

where  the  matrices  A  (iJ.^  (t))  and  B  (p^  (r))  have  the  form 

_  3  _  . 

A  (p^  (0)  =  A  0  +  X  A,-  (p ' (t )),  (2.3.3b) 

i=i 

B(iL,it))  =  Bo  +  B^(\i^{t)),  (2.3.3c) 

where  the  matrices  A,  (-),  /  =  1,2,3,  and  5  3(0  are  nonlinear,  continuously  differentiable 
functions  of  the  damping  coefficients  Pc(r),  and  Aq  and  Bq  defined  as  in  Section  2.2. 

The  matrices  A^(p^(f)),  /=1,2,3,  have  the  property  that  A,(0)  =  0  and 
Aj(Pc(t)) Aj  as  Pp(r)->oo.  Also,  53(0)  =  0  and  B3(Pc(r)) ->  B3  as  p^(r)->oo. 
Hence,  in  the  limit,  as  the  damping  coefficients  corresponding  to  the  clamps  go  to  0  or 
the  model  given  in  (2.3.3a)  becomes  identical  to  that  of  (2.2.4a)  with  the  corresponding 
clamp  status  set  to  0  or  1  respectively. 

Since  the  damping  vector  p^  (c )  is  continuous,  the  matrices  A  (p^,  (•))  and  B  (p^  (•)) 
are  continuous  functions  of  time,  and  by  the  fundamental  theorem  of  differential  equa¬ 
tions,  we  see  that  both  the  state  trajectories  z  (r )  and  their  derivatives  z  (r )  are  continuous 
for  all  t  G  [0,Tf]. 

Two  issues  still  remain  before  the  new  model  is  completely  specified.  These  are 
choosing  an  appropriate  clamp  transition  time,  T^,  and  an  exponent  r  for  the  damper 
response  law  in  (2.3.1a,b).  We  carried  out  numerical  experiments  using  the  model  in  this 
section  so  as  to  establish  reasonable  values  for  these  parameters.  Once  again,  we  used  a 
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sampling  time  T  =  0.2  seconds,  a  final  time  =  5  seconds,  the  earthquake  ground 
motion  supplied  by  the  El  Centro  1940  record,  and  a  command  signal  constructed  using  a 
simple  moving  horizon  control  law. 

exponent  for  damper  response  law  r.  We  set  =0.97,  which  gives  a 
sufficiently  long  clamp  transition  time  to  investigate  the  effects  of  different  values  of  r  in 
(2.3. 1  a, b).  We  found  that  r  =  4  gave  the  most  gradual  transition,  with  smaller  values  of 
r  resulting  in  most  of  the  smoothing  action  occurring  in  the  first  fifth  of  the  transition 
time,  while  larger  values  resulted  in  the  smoothing  action  occurring  in  the  last  fifth  of  the 
transition  time. 

transition  time  .  To  determine  if  the  value  7^  =  0.97  should  be  replaced  by 
another  value,  we  carried  out  simulations  using  7^  =  a7,  where  a  e  [0, 1),  for  several 
values  of  a.  It  was  observed  that  as  a  becomes  larger,  the  acceleration  peaks  in  the 
response  of  the  controlled  structure  become  smaller.  For  example,  see  Figures  2.3.1  - 
2.3.3,  where  Figure  2.3.1  shows  the  resulting  accelerations  using  a  =  0.0,  Figure  2.3.2 
shows  the  accelerations  for  a  =  0.5,  and  Figure  2.3.3  shows  the  accelerations  for  a  =  0.9. 

We  define  the  peak  acceleration  magnitude  vector  by 

Xpeak^i  max  max  lx,^(r)l,  max  f2  3  4al 

^  /  e  [0,7/1  <  e  [0,7/1  <  e  [0,7/1 

the  average  acceleration  vector  by 

\x,\t,c)\dt,j^'  \x,ht,c)\dt.  \x,\t,c)\dt),  (2.3.4b) 
and  the  maximum  interstory  drift  vector  x^  by 

Xj^i  max  IAx,^(r,c)l,  max  lAxAr,c)l,  max  I  Ax/(r,c)  I) .  o 

/  €  [0,7/1  t  €  [0,7/1  t  e  [0,7/1 

Referring  to  Table  2.3.1,  we  see  that  both  the  peak  and  average  accelerations  become 
smaller  as  a  increases. 


a 

1^91 

mm 

0.0 

575 

1025 

1200 

103.9 

164.6 

232.7 

0.5 

550 

800 

550 

113.0 

149.7 

128.4 

0.9 

440 

460 

410 

102.0 

109.5 

132.6 

Table  2.3.1.  Effect  of  a  on  Peak  and  Average  Accelerations. 

Hence,  our  results  seem  to  indicate  that  better  performance  in  terms  of  reducing  the 
accelerations  of  the  floors  of  the  structure  can  be  obtained  through  the  use  of  longer  tran¬ 
sition  times.  This  result  is  quite  important  since  it  indicates  that  it  may  not  be  necessary 
to  develop  faster  switching  dampers  in  order  to  improve  the  response  of  the  controlled 
structure.  In  the  remaining  simulations,  we  used  the  value  a  =  0.9. 
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In  order  to  determine  whether  the  simplified  model  considered  in  this  section  can 
always  be  used  in  numerical  experiments  in  place  of  the  original  system  model,  simula¬ 
tions  were  performed  using  both  models  and  identical  command  signals  and  ground 
motion  records.  The  command  signals  were  generated  using  a  simple  moving  horizon 
sampled-data  control  law  with  sampling  period  T  =  0.2  seconds,  and  the  earthquake 
ground  motions  are  given  by  the  El  Centro  (1940)  (EC),  Kern  County  (1952)  (KC),  Puget 
Sound  (1965)  (PS),  San  Fernando  (1971)  (SFO,  and  Western  Washington  (1949)  (WW) 
records.  A  comparison  of  the  peak  interstory  drifts  is  given  in  Table  2.3.2;  a  comparison 
of  the  peak  accelerations  is  given  in  Table  2.3.3;  and  a  comparison  of  average  accelera¬ 
tions  is  given  in  Table  2.3.4. 


simplified  model 

original  model 

Quake 

r  1 

y2 

r3 

^d 

r  1 
^d 

r2 

EC 

1.15 

1.90 

1.80 

1.13 

1.99 

1.89 

KC 

1.70 

2.90 

3.00 

1.79 

2.92 

3.22 

PS 

1.40 

1.65 

1.20 

1.42 

1.81 

1.42 

SF 

1.00. 

1.65 

1.00 

0.97 

1.73 

1.37 

1.40 

2.00 

1.40 

1.38 

2.13 

1.41 

Table  2.3.2.  Comparison  of  Peak  Interstory  Drifts  (cm ). 


simplified  model 

original  model 

Quake 

^peak 

■^peak 

^peak 

r'l 

■^peak 

r2 

^peak 

■*-peak 

EC 

625 

475 

525 

566 

358 

350 

KC 

525 

450 

475 

481 

396 

271 

PS 

390 

360 

375 

384 

361 

362 

SF 

380 

510 

400 

397 

380 

372 

WW 

275 

220 

380 

289 

243 

284 

Table  2.3.3.  Comparison  of  Peak  Accelerations  (cm  Isec^) 
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simplified  model 

original  model 

Quake 

..  1 
^avg 

^avg 

^avg 

^avg 

^avg 

^avg 

EC 

83.7 

79.6 

88.0 

82.4 

70.0 

82.0 

KC 

79.0 

74.1 

55.3 

67.1 

57.8 

44.5 

PS 

67.7 

73.3 

76.4 

72.6 

63.7 

79.8 

SF 

10.1 

62.3 

64.0 

65.4 

53.3 

54.9 

WW 

36.6 

37.8 

47.5 

34.8 

36.4 

45.3 

Table  2.3.4.  Comparison  of  Average  Accelerations  (cmlsec^) 

In  general,  the  simulations  using  the  simplified  model  exhibit  higher  accelerations 
and  lower  interstory  drifts  than  the  simulations  using  the  original  model.  These  differ¬ 
ences  are  too  significant  to  ignore,  and  hence  the  simplified  nKxiel  can  be  used  in  prelim¬ 
inary  experiments,  but  for  final  results,  the  original  model  must  be  used. 

3.  Estimation  of  Performance  Bounds 

Since  the  control  system  for  an  actively  controlled  variable  structure  cannot  impart 
energy  to  the  stracture,  it  is  not  clear,  a  priori,  what  is  its  potential  for  earthquake  damage 
mitigation.  Hence,  before  attempting  to  develop  control  laws,  we  will  establish  upper 
and  lower  bounds  on  its  performance.  The  upper  bounds  will  be  established  using 
optimal  control  laws  under  the  assumption  that  the  earthquake  ground  motion  is  known 
in  advance.  A  lower  bound  on  acceptable  performance  will  be  obtained  by  carrying  out  a 
minimax  optimal  design  that  assigns  fixed  values  to  the  variable  dampers,  on  the  basis  of 
the  histories  of  a  fixed  number  of  earthquakes.  The  first  step  in  obtaining  perfomoance 
bounds  is  the  development  of  a  performance  criterion  that  captures  the  essential  perfor¬ 
mance  requirements  of  interstory  drift  and  acceleration  limitation,  and  that  is  natural  for 
use  in  moving  horizon,  as  well  as  instantaneous  optimal,  feedback  control  laws. 

3.1.  Development  of  a  Performance  Criterion 

We  will  assume  that  minimal  performance  requirements  for  a  structure  subjected  to 
earthquakes,  scaled  to  have  the  same  peak  velocity  as  the  El  Centro  (1940  SOQE)  earth¬ 
quake,  are  as  specified  in  Table  3.1.1,  which  states  maximum  allowable  deformations  and 
shear  forces  that  the  structure  can  withstand  before  it  breaks.  In  Table  3.1.1, 
denotes  the  maximum  allowable  interstory  drift,  in  cm,  denotes  the  maximum 

allowable  shear  force,  in  N,  and  denotes  the  maximum  allowable  force  on  the 

damper,  in  N. 
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Story  1 

Story  2 

Story  3 

^max 

x}  <  2.02 

x/  <  2.46 

xj  <  1.89 

8.33  X  10^ 

7.448  X  10^ 

5.782  X  10^ 

^^max 

1.372  X  10^ 

1.372  X  10^ 

1.372  X  10® 

Table  3.1.1.  Design  Performance  Specifications 

In  addition,  the  accelerations  of  the  floor  masses  must  be  kept  as  small  as  feasible,  so  as 
to  help  protect  the  contents  of  the  structure  during  an  earthquake. 

The  above  requirements  are  best  captured  by  a  minimax  type  criterion.  However, 
there  are  no  commonly  used  feedback  laws  based  on  a  minimax  criterion,  and,  in  addi¬ 
tion,  it  would  be  rather  cumbersome  to  implement  a  moving  horizon  control  law  based 
on  a  minimax  criterion.  Hence  we  will  use  a  weighted  energy  criterion  that  reflects  these 
requirements  indirectly,  but  which  relates  naturally  to  common  feedback  laws.  The 
weights  will  be  tuned  to  make  the  energy  criterion  represent  the  performance  require¬ 
ments  as  well  as  possible. 

Given  a  ground  motion  vector  e  (t),  the  fact  that  the  state  of  the  variable  structure  is 
at  a  time  ^  0,  and  a  damper  coefficient  modulation  function  |j,(f),  we  will  denote 
the  resulting  response  of  the  structure,  determined  either  by  (2.1.1a,b)  or  by  (2.3.2a,b),  by 
z  (r ;  z^ ,  ,  p,  e ).  With  each  story  of  the  structure,  at  time  r  >  ,  we  will  associate  four 

energy  functions,  all  of  which  depend  on  the  variable  damper  modulation  functions  and 
the  ground  motion.  Let  be  the  kinetic  energy  in  i—th  story  of  the 

main  structure,  El  p  {t ;  z^ ,  ,  p,  c ),  the  potential  energy  in  i  -th  story  of  the  main  struc¬ 

ture,  El jy (t ; z^,T^,\i,e ),  the  potential  energy  in  the  i -th  story  of  the  auxiliary  struc¬ 
ture,  and  ),  the  kinetic  energy  in  the  i-th  story  of  the  auxiliary  struc¬ 

ture.  Then,  irrespective  of  whether  model  (2.1.1a,b)  or  model  (2.3.2a,b)  is  used,  we  see 
that  these  are  given  by 


;Zo.Toy\y.e)^VzM'lxi{tf,  i  =  1,2,3, 

(3.1.3a) 

Eijyit  ;z^,T^,\i,e)^V2Ki(xi-xi^'^f,  i  =  1,2,3, 

(3.1.3b) 

i  =  1,2,3, 

(3.1.3c) 

Eljpit  ;z„T„[i.e)kVzKi{x'l-xi^^f,  /  =1,2,3. 

(3.1.3d) 

We  will  use  as  our  performance  criterion  a  weighted  energy  function  of  the  form 

3  3 

E(t  ;z^,To,\L,e)^  E  +  E  ^ij>Elj,it  ;z^,T^,\L,e),  (3,i.3e) 

i=i  i=i 

where  the  weights  >0,i  =  1, 2, 3. 
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The  parameters  w  W2j> ,  ,  w  W2jc>  ^3,k)  enable  us  to  fine-tune  the 

response  of  the  controlled  structure.  Thus,  increasing  the  potential  energy  parameters 
Wij,  tends  to  decrease  corresponding  interstory  drifts,  while  increasing  the  kinetic  energy 
parameters  tends  to  reduce  the  velocity  of  the  i-th  story  mass. 

Using  the  cruder  model  (2.3.2a,b),  we  performed  a  set  of  numerical  experiments  to 
determine  the  weights  for  the  energy  function  in  (3.1.3e)  that  result  in  good  performance 
of  the  controlled  structure.  This  involved  performing  simulations,  using  Moving  Horizon 
Control  Law  4.1,  of  the  response  of  the  controlled  structure  to  the  the  1940  El  Centro 
(comp  SOOE)  earthquake  record  for  various  weights.  For  these  experiments,  the  initial 
time  was  =  0,  the  initial  state  =  0,  the  sampling  time  T  =  0.2  seconds,  and  the  final 
time  Tf  =  5  seconds. 

Numerical  experiments  were  performed  for  several  values  of  weighting  factors.  A 
comparison  of  the  resulting  peak  interstoiy  drift  values  for  each  of  the  weighting  factors 
is  given  in  Table  3.1.2. 


weighting  factor  w 

interstory  drifts  (cm) 

r  1 

r2 

r3 

Xd 

(1 , 1 . 1 , 1 , 1 , 1) 

1.85 

1.40 

4.20 

(3, 3, 3, 1,1,1) 

1.20 

2.20 

2.10 

(3 , 3.5 , 3.5 , 1 , 1) 

1.30 

1.80 

1.95 

(3, 2.5, 4, 1,1) 

1.15 

1.90 

1.80 

Table  3.1.2.  Peak  Interstory  Drifts  vs.  Weighting  Factor. 

From  these  results,  we  see  that  changing  the  weights  on  the  potential  energy  terms 
directly  affects  the  resulting  interstory  drifts  on  the  conesponding  floors.  The  use  of 
w  =  (3, 2.5, 4, 1, 1)  enabled  the  response  of  the  structure  to  meet  the  interstory  drift  con¬ 
straints  given  in  Table  3.1.1,  but  the  use  of  larger  weights  on  the  potential  energy  terms 
than  those  shown  in  the  above  results  did  not  enable  us  to  achieve  better  results.  This 
seems  to  indicate  that  the  amount  of  disturbance  rejection  possible  with  a  variable  struc¬ 
ture  system  is  limited. 

The  effectiveness  of  the  weights  w  =  (3, 2.5, 4, 1, 1),  was  verified  by  performing 
numerical  experiments  on  four  additional  earthquake  records.  In  addition  to  the  El  Cen¬ 
tro  record,  we  used  the  Kem  County  (Taft)  1952  (comp  S69E)  record,  the  Puget  Sound 
1965  (comp  S86W)  record,  the  San  Fernando  (Pacoima)  1971  (comp  S74W)  record,  and 
the  Western  Washington  1949  (comp  N86E)  record.  Since  these  records  represent  earth¬ 
quakes  of  varying  intensities,  the  ground  motion  records  needed  to  be  scaled.  It  is  com¬ 
mon  practice  to  scale  the  records  so  that  the  peak  velocities  are  equal  for  all  of  the 
records.  For  bur  numerical  experiments  throughout  this  report,  we  used  the  El  Centro 
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record  as  a  standard  and  scaled  the  other  earthquakes  so  that  their  peak  velocities  were 
equal  to  that  of  the  El  Centro  record.  Since  these  experiments  were  being  run  for  only  5 
seconds,  only  the  first  5  seconds  of  each  record  were  considered  when  determining  the 
scaling  factors.  This  resulted  in  the  scaling  factors  shown  in  Table  3.1.3. 


Earthquake 

Scaling  Factor 

KC 

1.825 

PS 

7.760 

SF 

0.563 

WW 

2.220 

Table  3.1.3.  Ground  Motion  Scale  Factors. 

The  scaled  ground  acceleration,  velocity,  and  position  for  each  of  the  earthquakes  are 
shown  in  Figures  3.1.1  -  3.1.15.  The  resulting  experiments  indicate  that  the  use  of 
w  =  (3, 2.5,4, 1, 1)  provides  reasonable  performance,  with  the  interstory  drift  constraints 
being  met  for  all  but  the  structure  subjected  to  the  Kem  County  earthquake. 

3.2.  Upper  Performance  Bounds:  Optimal  Controls 

From  now  on  we  will  use  only  the  system  model  (2.1.1a,b).  The  five  optimal  con¬ 
trol  problems  that  we  will  solve  differ  only  by  the  ground  motion  and  ground  velocity, 
obtained  from  an  earthquake  record,  appearing  in  the  system  model  (2.1.1a,b). 

the  optimal  control  problem.  Since  we  will  also  need  to  use  the  optimal  control 
algorithm  in  a  moving  horizon  control  law,  we  begin  by  stating  a  more  general  optimal 
control  problem  than  is  needed  for  this  section  alone.  Let  0  <  7^  <  Ty  <  oo  be  the  initial 
and  final  times.  Let  z  (r ;  ,  7^ ,  |X,  e )  denote  the  response  of  the  structure  model 

(2.1.1a,b),  with  initial  state  z^  at  time  7^  to  a  ground  motion  e  (f ),  when  controlled  by  the 
damper  modulation  function  {i(r ),  let 

M=  {peLl[7^,7^]  I0<4‘(r)<n„33„  /  =  1,2, 3,  f  e  [0,7^] } ,  (3.2.1a) 

and  let  <!>.•  M  ]R  be  defined  by 

d)(|i)  4  L E (r  ;z^,T^,[L,e)dt,  (3.2.1b) 

i  Q 

where  E{t ;  z^ ,  7^ ,  |X,  e )  is  defined  as  in  (3. 1 .3e).  To  compute  an  optimal  damper  modu¬ 
lation  function  (control)  p.(r),  we  propose  to  solve  the  following  optimal  control  prob¬ 
lem: 


3.2.  Upper  Performance  Bounds:  Optimal  Controls 


In  order  to  obtain  a  numerical  solution  to  P,  this  problem  must  be  replaced  by  a  con¬ 
sistent  finite  dimensional  approximation  [Pol.4].  To  obtain  such  an  approximation,  we 
replace  the  time  interval  by  a  grid  of  N  equally  spaced  time  points  =  k A, 

k  =  0, 1,...,  N,  where  A  =  [Tf  and  assume  that  the  the  controls  are  constant 

between  the  grid  points.  Assuming  without  loss  of  generality  that  =  0,  results  in  con¬ 
trols  having  the  finite  parametric  representation 


whereM(^)G  IR^.for^  =  l,2,..,N,and 


(3.2.2a) 


7t(f)  ^ 


1,  for  0  <  f  <  A, 
0,  otherwise . 


(3.2.2b) 


Hence  the  sequence  {u(k)]j^=i  completely  determines  the  control  ii(f).  Tlte  set  of 
admissible  discrete  control  sequences  u  =  {u(k)]^-i,  that  correspond  to  controls  in  the 
set  M,  is  given  by 


{«  =  {«(it)}f=i  lM(A:)e  V,  k  = 


(3.2.2c) 


where 

V4  {vGR3|o<v^<m^3^,  y  =  1,2,3).  (3.2.2d) 

At  the  same  time,  the  system  equations  (2.1.2)  are  replaced  by  the  backwards  Euler 
integration  formula: 

3  . 

F((;!:+1)A)  =  I(kA)  +  AiA  +  X  M*((ik-l-l)A)5,)F((it+l)A)  +  ADe((it+l)A),(3.2.3a) 

i=l 


where  due  the  the  linearity  of  the  system  with  respect  to  the  state  variables,  the  implicit 
difference  equations  can  be  solved  to  yield  the  explicit  difference  equation 

z((;k+l)A)  =  (/  -  A(A  +  £  M‘((A:-H)A)fi,))-i(F(^A)  + ADe((A:+l)A) .  (3.2.3b) 

i=l 

The  backwards  Euler  integration  method,  which  is  the  simplest  implicit  numerical 
integration  method,  is  used  here  since  for  large  values  of  |J.(t ),  the  state  equations  (2.1.2) 
become  stiff. 


Let  / ;  U  — »  R  be  defined  by 

f(u)^AJ^E(kA;Zo,T^,u,e),  (3.2.4) 

ifc=i 

where  E  defined  as  in  (3.1.3e),  but  with  ¥(kA)  replacing  z(kA)  in  equations  (3.1.3a-d). 
The  finite  dimensional  approximating  optimization  problem  is  then  given  by 

p„  tnin/(«).  (3.2.5) 

The  gradient  of  the  cost  function,  V/  (u )  is  an  element  of  R^^ .  To  retain  an  exact 
correspondence  to  the  scalar  product  and  norm  in  L2,  we  define  the  scalar  product 
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^Rby 

N 

ix,y)M  =^ZXi,yi,  (3.2.6a) 

i=l 

where  x  =  { jc,- } /Ij,  y=  and  g  R^.  The  corresponding  norm 

l  lyy ;  R^^  — >  R  is  then  given  by 

lhl^=  {h,h)isf  .  (3.2.6b) 

the  optimal  control  algorithm.  The  finite  dimensional  optimization  problem 
can  be  solved  numerically  using  a  number  of  optimization  algorithms.  We  will  use  the 
following  generalization  of  the  Armijo  gradient  method,  which  we  have  developed  for 
optimization  problems  with  simple,  convex  constraints,  as  in  (3.2.5)  (see  [Bak.l,  Pol.l, 
Pol.2]). 

Constrained  Armijo-Gradient  Algorithm  3.2.1. 

Parameters.  a,  P  g  (0, 1),  y>  0,  n,-  =  number  of  iterations. 

Data.  T^,Tf,z^e  R^^  e(.)^  =.  { u^ikA) )  g  U. 

Step  0.  Set  /  =  0. 

Step  1.  Compute  the  search  direction: 

hi  =  { {hi  (A),  hi  (2A), . . hi  (NA)}  4  arg  min  {  ( V/  (m,-  ),h)i^-\-  ^Ih  1^  I 

u{kA)  +  h{kA)s  V,  j  =  1,2,3,/:  =  l,2,...,iV  }.  (3.2.7a) 

•  Step  2.  Compute  the  step  size: 

Xi  =  max  \f{ui  +p'/i(M,))-/(M,)<ap^(V/(M.),/i(M,)V, 

I  e  IN 

uikA)  +  ^‘h{kA)e  V,  j  =  1,2,3,  k  =  l,2,...,iV  ).  (3.2.7b) 

Step  3.  Set  Ui+i  =  «,•  +  X,-  hi ,  replace  i  by  i +1 ,  and  if  i  <  rii ,  go  to  Step  1.  □ 

The  parameter  y  is  related  to  the  second  derivative  of  the  cost  function  in  the  fol¬ 
lowing  way: 

inf  ol/„„  (m  )]  <  Y  <  sup  aI/„„  (u )],  (3.2.8a) 

where  (5[Q],  c[Q]  denote  the  smallest  and  largest  singular  values,  respectively,  of  the 
matrix  Q .  Since  the  second  derivative  matrix  (u )  is  quite  tedious  to  compute,  and 
since  the  value  of  y  only  has  to  be  in  the  range  between  the  maximum  and  minimum 
singular  values,  an  acceptable  value  may  be  computed  simply  setting 
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for  any  m  e  U  and  8m  e  such  that  m  +  5m  €  U. 


(3.2.8b) 


One  final  note  on  Algorithm  3.2.1  is  that  the  stopping  criterion  may  be  changed  to 
ensure  that  the  algorithm  continues  until  it  approaches  a  local  minimum  rather  than  just 
for  a  specified  number  of  iterations.  For  example,  given  some  e  >  0,  temainate  the  itera¬ 
tion  when 

-{Vf(Ui),hi)<e.  (3.2.9) 

In  other  words,  when  the  directional  derivative  in  the  direction  of  the  search  direction 
becomes  small  (in  magnitude),  additional  iterations  will  not  produce  any  significant 
decreases  in  the  value  of  the  cost  function,  so  the  iteration  should  be  terminated. 


Initialization  algorithm  Optimal  control  problems  require  considerable  comput¬ 
ing  time  to  solve,  with  the  solution  time  depending  strongly  on  the  choice  of  the  initial 
control.  We  therefore  propose  to  solve  the  optimal  control  problems  in  two  stages.  The 
first  consists  of  solving  a  simplified  problem,  in  which  the  control  is  required  to  be  con¬ 
stant  for  all  times,  for  an  initial  control  to  be  used  in  the  second  stage  by  Algorithm  3.2.1. 
Referring  to  (3.2.5),  we  see  that  the  simplified  problem  must  have  the  form 

min  {/ (m)  I  Ml  =  M2  =  . . .  =  Mjv  } ,  (3.2.10a) 

«  e  U  '  ' 

or,  equivalently, 

po  min/(v),  (3.2.10b) 

where 

/(v)4/(M(v)).  (3.2.10c) 

with  u(y)=  { V  )  j^i,  i.e.,  Mjt(v)  =  v  for  all  k. 

Given  the  form  of  /  (v ),  we  see  that  we  can  compute  its  gradient  by  simply  using 
the  chain  rule.  The  result  is 

V/(v)  =  M,(v)^V/(M(v)),  (3.2.11) 

where  u^(v)^  =  [73,73, ...,73]  is  a  3x3N  matrix,  with  1 2  the  3x3  identity  matrix. 
Hence,  when  applied  to  problem  ,  Algorithm  3.2. 1  assumes  the  form 

Constrained  Armijo-Gradient  Algorithm  3.2.1’. 

Parameters.  yvij>,W2j>, ^3, if  > 0,  a, P  e  (0, 1), y > 0,  «,•  =  number 
of  iterations. 

Data.  VqS  V. 

Step  0.  Set  j  =  0. 


-20- 


3.2.  Upper  Performance  Bounds:  Optimal  Controls 


Step  1.  Compute  the  search  direction: 

hi  4  arg  min  {  ( V/ (vj ),  h  +  -Jll/i  Ijv  } .  (3.2. 12a) 

h  eV -Vi  ^ 

Step  2.  Compute  the  step  size: 

Xi  =  max  {  p*  l/(v,  +p^/j,))-/(v,)<aP'(V/(u,),/t(M,)V,,  v,  +p^/i,-  e  V}  .(3,2.12b) 

/  €  N 

Step  3.  Set  Ui+i  =  Ui  +  X,-  hi ,  replace  i  by  i  +1 ,  and  if  i  ^  n,- ,  go  to  Step  1 .  □ 

Problem  Pjy  takes  less  time  to  solve  than  problem  P;^,  because  the  dimension  of  the 
search  direction  calculation  in  Step  1  of  Algorithm  3.2.1’  is  N  times  smaller  than  in 
Algorithm  3.2.1.  In  addition,  our  computational  experience  indicates  that,  in  general, 
fewer  iterations  are  necessary  to  solve  problem  P^^  than  problem  . 

computational  experiments  The  parameters  used  in  our  optimal  control  experi¬ 
ments,  both  for  Algorithm  3.2.1  and  Algorithm  3.2.1’  are  given  in  Table  3.2.1. 


w 

A  (sec ) 

Tf  (sec) 

a 

p 

Y 

(3, 2.5, 4, 1,1,1) 

0.004 

5 

0.5 

0.75 

1 X 10"^® 

Table  3.2.1.  Parameters  Used  in  Optimal  Control  Experiments 

•  initialization  Optimal  constant  controls  (i.e.,  fixed  values  for  the  damping 
coefficients),  for  use  as  initial  controls  by  Algorithm  3.2.1,  were  computed  using  Algo¬ 
rithm  3.2.  r  for  each  of  the  five  scaled  earthquake  ground  motion  records.  The  parameter 
Y  was  selected  by  means  of  an  expression  analogous  to  (3.2.8b).  The  fact  that  gamma  is 
quite  small  indicates  that  the  cost  function  /  (v )  is  close  to  being  linear.  All  of  the  com¬ 
putations  with  Algorithm  3.2.1’  required  no  more  than  10  iterations  for  the  termination 
test  I  ( V/  (v,  ),  hi )  / 1  <  0.01/  (vq)  to  be  satisfied. 

The  resulting  optimal  fixed  damping  coefficient  values  are  shown  in  Table  3.2.2. 
We  see  that  the  optimal  fixed  damping  coefficient  for  the  third  floor  is  largest  and  that  the 
optimal  fixed  damping  coefficients  for  the  first  and  second  floors  are  quite  close.  This 
intuitively  makes  sense  since  the  weight  on  the  potential  energy  terms  for  the  cost  func¬ 
tion  is  largest  for  the  third  floor,  while  weights  the  first  and  second  floors  are  smaller,  and 
almost  equal. 
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Earthquake 

El  Centro 

Kern  County 

7.9  X  10^ 

7.6  X  10^ 

1.80  X  10^ 

Puget  Sound 

9.1  X  10^ 

8.7  X  10^ 

1.71  X  10'^ 

San  Fernando 

1.72  X  10^ 

1.63  X  10^ 

3.15x10^ 

Western  Washington 

1.43  X  10“^ 

1.44x10^ 

3.20x10”^ 

Table  3.2.2.  Initial  Design:  Optimal  Fixed  Damper  Values 

•  optimal  controls.  For  Algorithm  3.2.1,  we  used  as  initial  controls  Uq,  the  constant 
controls  computed  above.  In  contrast  with  our  initialization  computations  using  Algo¬ 
rithm  3.2.  r,  the  computation  of  optimal  controls  using  Algorithm  3.2.1  required  as  many 
as  25  iterations  to  reach  the  termination  test.  Also,  as  expected,  each  iteration  of  Algo¬ 
rithm  3.2.1  took  somewhat  longer  than  each  iteration  of  Algorithm  3.2.1’.  Our  results 
indicate  that  in  general,  when  the  earthquake  ground  motions  are  small,  the  optimally 
modulated  damping  coefficients  remain  close  to  the  initial  fixed  values,  but  when  the 
ground  motion  becomes  more  severe,  the  optimally  modulated  damping  coefficients  can 
deviate  significantly  firom  the  initial  fixed  values.  We  see  that  for  each  earthquake,  the 
optimally  controlled  structure  easily  meets  the  interstory  drift  specifications  in  Table 
3.1.1.  A  summary  of  the  results  of  each  of  the  computations  will  be  given  in  Section  5 
along  with  a  detailed  discussion  of  these  results.  In  the  interest  of  keeping  the  size  of  this 
report  reasonable,  we  have  included  a  full  set  of  plots  for  only  the  computations  which 
use  the  El  Centro  ground  motion  record.  The  plots  are  given  in  Figures  3.2.1 -3.2.6, 
where  Fig.  3.2.1  shows  the  interstory  drifts.  Fig.  3.2.2  shows  the  floor  accelerations.  Fig. 
3.2.3  shows  the  modulated  damping  coefficients.  Fig.  3.2.4  shows  the  energy  dissipation 
in  the  variable  dampers.  Fig.  3.2.5  shows  the  shear  force  on  the  frame  of  the  structure, 
and  Fig.  3.2.6  shows  the  force  on  the  variable  dampers. 

3.3.  Lower  Performance  Bounds:  Minimax-Optimal  Fixed  Structure 

To  justify  the  additional  cost,  a  controlled  variable  structure  must  perform  better 
than  any  well  designed,  comparable  fixed  structure.  We  propose  to  obtain  minimax- 
optimal  fixed  values  for  the  damping  coefficients  with  respect  to  the  Kern  County,  Puget 
Sound,  San  Fernando,  and  Western  Washington  earthquakes.  The  El  Centro  earthquake 
will  not  be  included  in  our  minimax  design,  so  that  we  can  later  compare  the  responses  of 
the  optimally  controlled  structure  and  minimax-optimal  designed  fixed  structure  on  an 
earthquake  that  was  not  included  in  the  design  of  the  fixed  structure. 

formulation  of  the  minimax-optimal  design  problem.  Again  the  continuous 
system  dynamics  must  be  di^retized  by  the  integration  formula,  as  described  in  Section 
3.2.  For  j  =  1,2, 3,4,  let  f^iv)  be  defined  as  in  (3.2.13c),  with  the  superscript  j 


-22- 


3.3.  Lower  Performance  Bounds:  Minimax-Optimal  Fixed  Structure 


denoting  the  particular  earthquake  exciting  the  dynamic  equations  (3.2.6a),  as  follows: 
7  =  1  corresponds  to  Kem  County,  j  =  2  corresponds  to  Puget  Sound,  7=3  corresponds 
to  San  Fernando,  and  7=4  corresponds  to  Western  Washington.  Then  we  propose  to 
obtain  fixed  values  for  the  dampers  by  solving  the  problem 


T\mm 

min  max  (/-^(v) ) , 

(3.3.1) 

V  E  V  7  €  m 

where m  =  {1,2, 3,4}.  Let 

\|/(v)  =  max  /^  (v) . 

7  6  m 

(3.3.2) 

The  problem  can  be  solved  by  means  of  the  following  constrained  minimax 
algorithm,  which  is  a  direct  extension  of  Algorithm  3.2.1’  (see  [Bak.l,  Pol.l,  Pol.2]). 

Constrained  Minimax  Algorithm  3.3.1. 

Parameters.  w  =  (wij>,W2j>,W2j>,Wi^,W2,K>^-i,K)y  a,Pe(0, 1),  y>0,  n,- = 
number  of  iterations. 

Data.  V  e  V. 

Step  0.  Set  I  =  0. 

Step  1.  Compute  the  search  direction  ft,-  and  value  of  the  optimality  function  0,- : 

hi  4  arg  min  max  { fHvi )  -  \jr(Vj )  +  ( V/-'  (v,-  ),h)+^lhf} ,  (3.3.3a) 

heV-Vijem  2 

0,  =  min  max  { p  (v,- )  -  \jr(v,- )  +  ( Vp  (v;  \h)+^lh  f  } .  (3.3.3b) 

€  V“V,  7  G  m  ^ 

Step  2.  Compute  the  step  size: 

Xi  =  max  {  I  \|/(v,  +  p* hi )  -  v(v, )  <  ap^  0(v,- ),  v,-  +  P^  /i,-  e  V } .  (3.3.3c) 

it  €  N 

Step  3.  Set  =  v,-  +  A,,-  hi ,  replace  i  by  i  +1,  and  if  i  <  rii ,  go  to  Step  1.  □ 

An  alternative  stopping  criterion  for  this  algorithm  is  to  stop  iterating  when 
1 0j  I  <  e,  where  e  >  0  is  small.  This  indicates  that  the  current  vector  of  damping 
coefficients,  v,- ,  is  close  to  a  local  minimizer  of  the  max  function  \j/(  ). 

minimax-optimal  design  of  fixed-structure.  We  set  a  =  0.9,  P  =  0.75,  and,  as 
in  the  previous  experiments,  y  =  1  x  10“^®.  The  initial  damping  coefficient  vector  was 
chosen  to  be  vq  =  (5  x  10^5  x  10^,5  x  10^),iV5/m.  Algorithm  3.3.1  required  only  7 
iterations  to  compute  a  v,-  such  that  the  stopping  criterion  I0(v,  )l  <0.01\|/(vo)  was 
satisfied.  The  resulting  damping  coefficient  values  were 

^  =  (0.90  X  10'^,0.87  X  lO"^,  1.70  x  lO"^)  Nsedm, 

which  happen  to  be  nearly  identical  to  the  fixed  optimal  values  computed  for  the  Puget 
Sound  earthquake,  in  the  initialization  stage  of  the  optimal  control  computations  of 


-23- 


3.3,  Lower  Performance  Bounds:  Minimax-Optimal  Fixed  Structure 


Section  3.2.  The  response  of  the  minimax  designed  structure  was  computed  for  each  of 
the  five  earthquakes  and  a  detailed  evaluation  will  be  carried  out  in  Section  5.  A  full  set 
of  plots  are  included  for  the  response  of  the  minimax  designed  structure  to  the  El  Centro 
quake.  Figure  3.3.1  shows  the  interstory  drifts.  Fig.  3.3.2  shows  the  floor  accelerations. 
Fig.  3.3.3  shows  the  energy  dissipation  in  the  optimally  chosen  dampers.  Fig.  3.3.4 
shows  the  frame  shear  forces,  and  Fig.  3.3.5  shows  the  force  on  the  dampers. 


4.  Moving  Horizon  Feedback  Controi  Laws 

We  will  now  describe  two  moving  horizon  control  laws  that  we  evaluated  with 
respect  to  our  performance  bounds.  The  first  is  a  sampled-data  relay  type  law  that  con¬ 
trols  the  dampers  in  on-off  fashion,  while  the  second  is  a  continuous  modulation  type 
law.  Both  of  these  moving  horizon  control  laws  must  be  considered  to  be  conceptual, 
since  they  requires  knowledge  of  the  ground  motion  one  horizon  ahead  of  the  present 
time,  which  may  not  be  available  in  a  real  time  control  system. 

sampled-data  moving  horizon  feedback  control  law  Our  sampled-data  moving 
horizon  control  law  is  based  on  a  sampling  period  T.  At  time  r  =  ^T,  k  e  IN  (i.e.,  at  the 
beginning  of  the  k-th  sampling  interval),  the  control  law  chooses  an  element  e  C 
that  determines,  via  (2.3.1a-d)  the  action  of  the  dampers  over  the  k-th  sampling  interval, 
so  as  to  minimize  the  cost  function 

E i{k-k\)T  ; z (kT), kT,\ic,e).  (4.1) 

For  i  =  1, 2, . . .,  8,  let  Ji  denote  the  elements  of  C.  Then  our  sampled-data  moving  hor¬ 
izon  control  law  may  be  formally  presented  as  follows: 

Moving  Horizon  Control  Law  4.1. 

Parameters.  w  =  (wij,.W2j>,W2jp,Wij^,W2,Kf^3,K)' 

Data.  Cqe  C. 

Step  0.  Set  k  =  0. 

Step  1.  At  t=kT,  compute  E((k+1)T  ;z(kT),kT,ii^  ,e),  j=l,2, ...,8,  where 

Cf 

diit)  =  Yi  for  r  e  [kT,(k+l)T)  and  ?,(r)  =  c(t)  for  r  e  [0,ytr). 

Step  2.  Set  Cf^  =  arg  min  E((k+1)T  ;z(kT),kT ,\i^  ,e),  and  for  t  e  [kT ,  (k+l)T),  set 

Yi  e  c  c . 

c(r)  =  c^. 

Step  3.  Replace  k  by  ^  +  1  and  go  to  Step  1.  □ 

With  the  sampling  time  set  at  T  =0.2  sec.,  the  weight  w  =  (3,2.5,4, 1, 1),  as 
selected  in  Section  3.1,  ,the  dynamics  determined  by  (2.1.1a),  and  the  damper  modulation 
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law  by  (2.3.1a-d),  we  simulated  the  response  of  the  structure,  controlled  by  Moving  Hor¬ 
izon  Control  Law  4.1,  with  Cq  =  (0,0,0),  subject  to  the  ground  motions  specified  by  the 
first  five  seconds  of  our  five  earthquake  records,  scaled  as  specified  Section  3.1.  Our 
results  indicate  that  the  the  response  of  the  controlled  structure  to  the  Puget  Sound,  San 
Fernando,  and  Western  Washington  earthquakes  meets  the  interstory  drift  constraints 
given  in  Table  3.1.1,  the  response  to  the  El  Centro  earthquake  slightly  violates  the  con¬ 
straints,  and  the  response  to  the  Kern  County  earthquake  significantly  violates  the  con¬ 
straints.  In  Section  5,  the  performance  of  the  control  law  for  each  of  the  five  earthquakes 
is  summarized  and  evaluated.  For  the  structure  subjected  to  the  the  El  Centro  earth¬ 
quake,  a  full  sets  of  plots  is  included  in  Figures  4.1  -  4.6,  where  Fig.  4.1  shows  the  inters- 
tory  drifts.  Fig.  4.2  shows  the  floor  accelerations.  Fig.  4.3  shows  the  variable  damping 
coefficients  generated  by  Control  Law  4.1,  Fig.  4.4  shows  the  energy  dissipated  by  the 
variable  dampers.  Fig.  4.5  shows  the  frame  shear  forces,  and  Fig.  4.6  shows  the  forces  on 
the  variable  dampers. 

continuous  time  moving  horizon  feedback  control  law  Our  continuous  time 
moving  horizon  control  law  is  based  on  a  horizon  of  T  seconds.  At  time  t  =  kT,k  e  N, 
the  control  law  uses  the  Optimal  Control  Algorithm  3.2. 1  to  determine  an  optimal  control 

P-;k(f).  which  is  then  used  over  the  time  interval  [kT ,(k+l)T). 

Moving  Horizon  Control  Law  4.2. 

Parameters.  w  = 

Step  0.  Set  ^  =  0. 

Step  1.  At  time  t  =kT,  set  up  data  for  Algorithm  3.2.1:  ZQ  =  z(kT),  -  kT, 
Ty  =  (A:+l)r,  and  Uq. 

Step  2.  Use  Algorithm  3.2.1  to  compute  an  optimal  modulation  function  p.^(r),  with 
t  &[kT,i)c+\)T). 

Step  3.  Apply  the  modulation  function  flj^(t),  for  t  e  [kT,  (k+l)T). 

Step  4.  Replace  k  by  k  +  I  and  go  to  Step  1.  □ 

Note  that  the  above  Control  Law  does  not  require  a  specific  rule  for  selecting  Mq  in 

Step  1  above.  For  our  experiments,  we  set  Mo(0  =  0  for  k  =0,  and  UQ(t)  =  fr_i(kT),  for 
k>0. 

For  our  numerical  simulations,  we  once  again  used  the  weight  w  =  (3, 2.5,4, 1, 1), 
the  dynamics  determined  by  (2.1.1a),  and  the  ground  motion  given  by  each  of  our  five 
scaled  ground  motion  records.  We  set  the  horizon  length  T  equal  to  0.2  seconds.  Our 
results  indicate  that  the  response  of  the  controlled  structure  to  each  of  the  earthquakes 
easily  satisfies  the  interstory  drift  constraints  given  in  Table  3.1.1.  A  complete  evalua¬ 
tion  of  the  performance  of  Control  Law  4.2  is  given  in  Section  5.  For  the  El  Centro 
earthquake,  we  have  included  a  full  set  of  plots  in  Figures  4.7  -  4.12,  where  Fig.  4.7 
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shows  the  interstory  drifts.  Fig.  4.8  shows  the  floor  accelerations,  Fig.  4.9  shows  the 
variable  damping  coefficients  generated  by  Control  Law  4.2,  Fig.  4.10  shows  the  energy 
dissipated  by  the  variable  dampers.  Fig.  4.11  shows  the  frame  shear  forces,  and  Fig.  4.12 
shows  the  forces  on  the  variable  dampers. 

5.  Evaluation  of  Numerical  Results 

We  will  now  discuss  the  results  of  our  numerical  experiments.  Because  of  the  need 
to  perform  a  very  large  number  of  numerical  experiments,  to  reduce  the  total  computa¬ 
tion  time,  we  performed  experiments  using  only  the  first  5  seconds  of  each  earthquake. 
In  addition,  to  establish  the  validity  of  our  conclusions  based  on  these  “short  time” 
experiments,  we  performed  one  set  of  experiments  using  the  ground  motion  supplied  by 
the  first  20  seconds  of  the  El  Centro  earthquake.  In  this  section,  we  will  the  refer  to  the 
optimal  designs  of  a  fixed  structure,  obtained  in  Section  3.2,  using  Algorithm  3.2.1’,  as 
initializations  for  the  optimal  control  computations,  as  optimal  fixed  designs.  To  gain 
more  insight  into  our  evaluation  of  the  controlled  and  minimax  designed  structures,  we 
have  computed  the  response  of  the  optimal  fixed  designed  structures  subject  to  the 
respective  earthquakes.  The  results  are  summarized  and  evaluated  in  this  section.  In 
addition,  in  Figures  5.1  -  5.5,  we  have  included  a  set  of  plots  of  the  response  of  the 
optimal  fixed  design  for  the  El  Centro  earthquake.  Specifically,  Fig.  5.1  shows  the  inters¬ 
tory  drifts.  Fig.  5.2  shows  the  floor  accelerations.  Fig.  5.3  shows  the  total  energy  dissipa¬ 
tion  in  the  optimally  chosen  dampers.  Fig.  5.4  shows  the  frame  shear  force,  and  Fig.  5.5 
shows  the  force  on  the  dampers. 

In  the  tables  below,  EC  denotes  the  El  Centro  1940  (SOOE)  earthquake,  KC  denotes 
the  Kem  County  (Taft)  1952  (comp  S69E)  earthquake,  PS  denotes  the  Puget  Sound  1965 
(comp  S86W)  earthquake,  SF  denotes  the  San  Fernando  (Pacoima)  1971  (comp  S74W) 
earthquake,  and  WW  denotes  the  Western  Washington  1949  (comp  N86E)  earthquake. 

optimal  control  versus  minImax  design.  We  will  compare  the  simulation 
results  obtained  using  optimal  control  of  the  variable  structure  with  those  obtained  using 
the  fixed,  minimax  designed  structure,  so  as  to  evaluate  the  available  performance  gap. 

•  optimal  cost  gap.  We  begin  with  the  gap  between  the  upper  performance  bounds 
resulting  from  using  optimal  control  with  perfect  information  about  the  earthquakes,  and 
a  lower  acceptance  bound  based  on  our  minimax  design  of  fixed  damping  coefficient 
values. 

First,  we  will  examine  the  primary  gap  Gp ,  for  each  of  our  five  earthquakes,  in 
terms  of  the  integral  of  our  weighted  energy  function  E,  defined  in  (3.1.3e).  We  define 
the  gap  by  the  fractional  relationship 
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g  (5.1, 

fm 

where  /  (m  ),  defined  in  (3.2.4),  is  the  cumulative  weighted  energy  cost  resulting  from  the 
use  of  an  optimal  control,  for  a  particular  earthquake,  and  /  defined  in  (3.2.10c), 

is  the  cumulative  weighted  energy  cost,  for  a  particular  earthquake,  associated  with  the 
minimax  designed  structure. 

Referring  to  Table  5.1,  we  see  that  the  biggest  gap,  of  61%,  occurs  for  the  El  Centro 
earthquake,  for  the  remaining  earthquakes,  the  gap  varies  from  27%,  for  the  Kern  County 
earthquake,  to  41%  for  the  San  Fernando  earthquake. 


Optimal 

Control 

Minimax 

Design 

Fixed 

Design 

Moving 

Horizon 

S-D  Moving 
Horizon 

Quake 

fm 

mm) 

f  fd) 

/  (^mh) 

/  (^sd) 

EC 

2.77 

4.46 

4.41 

3.26 

4.93 

KC 

1.65 

2.10 

2.09 

1.83 

4.92 

PS 

1.93 

2.51 

2.51 

2.65 

3.56 

SF 

1.58 

2.22 

2.01 

2.11 

2.89 

WW 

1.21 

1.57 

1.49 

1.27 

2.08 

Table  5.1.  Comparison  of  Optimal  Cost  Function  Values  (Nms  x  10*) 


Next  we  will  examine  how  the  primaiy  gaps  translate  themselves  into  gaps  between 
the  important  response  parameters,  viz.,  peak  and  average  accelerations  of  the  floors, 
peak  interstory  drifts,  total  energy  dissipation  in  the  variable  dampers,  peak  forces  on  the 
variable  dampers,  and  peak  shear  forces  on  the  main  structure. 

•  peak  interstory  drifts.  We  see  from  Table  5.2,  that  the  peak  interstory  drifts  reflect  our 
choice  of  weights  on  the  potential  energy  terms  in  the  energy  function,  i.e.,  the  drifts  are 
highest  on  the  second  floor,  where  the  weight  on  the  potential  energy  was  the  smallest. 
For  the  El  Centro  earthquake,  the  peak  interstory  drift  of  the  minimax  design  is  64% 
larger  than  that  of  the  optimally  controlled  structure,  and  so  we  see  that  the  use  of 
optimal  control  produces  a  considerable  reduction  in  maximum  interstory  drift  over  the 
minimax  design  for  this  particular  earthquake.  However,  for  the  other  earthquakes,  the 
superiority  of  the  optimal  control  law  was  not  clear. 
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Optimal  Control 

Minimax  Design 

Optimal  Design 

Quake 

r  1 

Xd 

r2 

Xd 

„3 

Xd 

Xd 

r2 

Xd 

Xd 

Xd 

r2 

Xd 

„3 

Xd 

EC 

0.80 

1.45 

1.10 

1.30 

1.90 

1.75 

1.20 

1.90 

1.60 

KC 

0.65 

1.45 

1.00 

0.75 

1.25 

1.10 

0.90 

1.30 

1.00 

PS 

0.75 

1.18 

0.70 

0.65 

0.95 

1.00 

0.60 

1.00 

1.00 

SF 

0.62 

1.12 

1.08 

0.85 

1.30 

1.15 

0.75 

1.15 

1.15 

WW 

0.75 

1.30 

0.80 

0.65 

1.05 

0.90 

0.60 

0.85 

0.80 

Table  5.2.  Comparison  of  Peak  Interstory  Drifts  (cm ) 

•  peak  accelerations.  Referring  to  Table  5.3,  we  see  that  for  the  El  Centro  earthquake, 
the  minimax  designed  structure  exhibits  30  -  70%  larger  accelerations  than  the  one  con¬ 
trolled  by  optimal  control,  with  the  biggest  difference  occurring  on  the  top  floor.  How¬ 
ever,  for  the  remaining  earthquakes,  the  superiority  of  the  optimally  controlled  structure 
is  again  not  clear. 


Optimal  Control 

Minimax  Design 

Optimal  Design 

Quake 

r  1 

■^■peak 

^peak 

■^peak 

r  ^ 
^peak 

■''■peak 

■''•peak 

r* 

■^peak 

•^peak 

r-3 

•^peak 

EC 

290 

260 

280 

500 

410 

360 

460 

375 

340 

KC 

320 

275 

290 

350 

280 

280 

340 

275 

280 

PS 

260 

220 

330 

235 

255 

320 

230 

250 

325 

SF 

310 

300 

360 

315 

280 

210 

350 

310 

225 

WW 

280 

260 

290 

265 

220 

225 

310 

260 

240 

Table  5.3.  Comparison  of  Peak  Accelerations  {cm/sec^) 


•  average  accelerations.  Referring  to  Table  5.4,  we  see  that  the  optimally  controlled 
structure  sustained  10-20%  higher  average  accelerations  on  the  ground  floor  than  the 
minimax  design,  but  the  minimax  design  sustained  up  to  65%  higher  average  accelera¬ 
tions  on  the  top  two  floors  than  the  optimally  controlled  structure,  with  the  largest  differ¬ 
ence  occurring  for  the  El  Centro  earthquake. 
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Optimal  Control 

Minimax  Design 

Optimal  Design 

Quake 

^avg 

JC^ 

•^avg 

•3 

^avg 

^avg 

••2 

Xavg 

^avg 

^avg 

jc-2 

•^avg 

■^avg 

EC 

59.6 

57.5 

79.8 

98.5 

80.2 

72.8 

91.6 

75.4 

72.1 

KC 

38.0 

36.5 

50.3 

50.4 

43.4 

44.1 

45.1 

42.6 

44.2 

PS 

43.8 

50.5 

77.9 

64.0 

58.6 

70.6 

64.0 

58.6 

70.7 

SF 

47.2 

59.8 

54.4 

47.6 

49.9 

62.7 

52.2 

50.5 

WW 

35.3 

45.3 

35.8 

34.2 

41.3 

41.8 

39.3 

43.7 

Table  5.4.  Comparison  of  Average  Accelerations  (cmlsec^) 

•  variable  clamper  energy  dissipation.  In  Table  5.5,  we  denote  the  total  energy  dissi¬ 
pated  in  the  variable  damper  connected  to  the  ith  floor  of  the  structure,  up  to  time  =  5 
seconds,  by  Ej-.  Referring  to  Table  5.5,  we  see  that  a  significant  advantage  of  using 
optimal  control,  over  a  minimax  design,  is  that  it  results  in  considerably  lower  energy 
dissipation  in  the  dampers.  We  see  that  the  dampers  of  the  minimax  design  must  dissi¬ 
pate  up  to  2.5  times  more  energy  than  the  dampers  of  the  optimally  controlled  structure. 


Optimal  Control 

Minimax  Design 

Optimal  Design 

Quake 

Q| 

Et 

E^ 

E^ 

E^ 

E^ 

EC 

1.20 

3.75 

3.10 

2.60 

6.20 

4.60 

2.60 

6.44 

4.53 

KC 

0.70 

2.10 

1.20 

1.00 

2.40 

1.80 

1.06 

2.51 

1.62 

PS 

0.61 

1.50 

1.40 

0.90 

2.1 

1.70 

0.93 

2.10 

1.66 

SF 

0.41 

1.52 

0.81 

1.00 

2.35 

1.80 

0.73 

1.72 

1.24 

WW 

0.33 

0.87 

0.50 

0.46 

1.12 

0.92 

0.42 

0.96 

0.67 

Table  5.5.  Comparison  of  Total  Energy  Dissipation  in  Dampers  (Nm  xlO^) 

•  main  structure  shear  force.  In  Table  5.6,  we  denote  the  peak  shear  force  at  the  i-th 
floor  of  the  main  structure  by  Fj.  We  see  from  Table  5.6  that  for  the  El  Centro  and  San 
Fernando  earthquakes,  the  minimax  structure  exhibits  4%  -  65%  higher  structural  shear 
forces  than  the  optimally  controlled  structure,  but  for  the  other  earthquakes,  the  advan¬ 
tage  of  the  use  of  optimal  control  over  a  minimax  design  is  not  as  clear. 
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Optimal  Control 

Minimax  Design 

Optimal  Design 

Quake 

H 

D 

B 

O 

B 

B 

B 

D 

EC 

1.7 

2.8 

2.6 

2.6 

3.8 

4.3 

2.5 

3.7 

4.0 

KC 

1.4 

2.4 

2.9 

1.7 

2.5 

2.6 

1.8 

2.6 

2.4 

PS 

1.5 

2.3 

1.7 

1.2 

1.9 

2.5 

1.2 

1.9 

2.5 

SF 

1.2 

2.2 

1.7 

2.6 

2.8 

1.5 

2.3 

2.9 

WW 

1.4 

2.5 

2.0 

1.3 

2.0 

2.2 

1.2 

1.7 

2.0 

Table  5.6.  Comparison  of  Peak  Shear  Force  (N  x  10^) 

•  peak  clamper  force.  The  peak  force  on  the  damper  connecting  the  i-th  level  of  the 
main  and  auxiliary  structures  is  denoted  in  Table  5.7  as  F|i.  Referring  to  Table  5.7,  we 
see  that  for  the  El  Centro  earthquake,  the  peak  force  on  the  variable  dampers  is  25%-50% 
larger  for  the  minimax  design  than  for  the  optimally  controlled  structure,  but  for  the 
remaining  earthquakes,  the  controlled  structure  exhibits  up  to  30%  higher  forces  on  the 
dampers  than  the  minimax  design. 


Optimal  Control 

Minimax  Design 

Optimal  Design 

Quake 

B 

B 

B 

B 

B 

B 

F  ^ 

B 

B 

EC 

5.3 

8.5 

9.6 

7.5 

11.5 

14.0 

6.9 

10.4 

12.8 

KC 

5.2 

9.0 

10.0 

1031 

8.3 

10.0 

8.0 

9.8 

PS 

4.6 

Bfll 

3.6 

5.2 

6.8 

3.6 

5.2 

6.8 

SF 

5.1 

9.1 

10.0 

4.9 

9.6 

noil 

8.8 

10.5 

WW 

5.2 

7.9 

iBBI 

4.3 

6.3 

7.5 

5.1 

7.5 

8.8 

Table  5.7.  Comparison  of  Peak  Damping  Force  (N  x  10^) 

•  conclusions.  As  expected,  on  all  earthquakes,  the  optimally  controlled  structure 
results  in  smaller  cost  function  values  than  the  minimax  design,  with  the  gap  being  larg¬ 
est  for  the  El  Centro  earthquake.  Our  results  also  show  that,  with  the  exception  of  the  El 
Centro  earthquake,  the  primary  performance  gap  Gp  does  not  clearly  translate  into  a 
corresponding  gap  in  structural  response  performance.  Overall,  the  performance  gap 
between  the  realistic,  minimax  designed  stmcture  and  the  conceptual,  optimally  con¬ 
trolled  structure  is  fairly  small,  and  thus  it  may  be  quite  difficult  to  obtain  a  feedback 
control  law  tiiat  does  significantly  better  than  the  minimax  design. 
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evaluation  of  minimax  design.  Our  minimax  design  was  based  on  a  very  small 
number  of  earthquakes.  We  will  now  examine  how  well  the  minimax  design  performs 
relative  to  the  unattainable  design  of  a  fixed  structure  specifically  reconfigured  for  each 
of  the  earthquakes  used  in  the  minimax  design,  as  well  as  relative  to  a  fixed  optimal 
design  for  the  El  Centro  earthquake,  which  was  not  used  in  the  minimax  design. 

•  optimal  cost  gap.  First,  we  will  consider  the  gap  in  the  optimal  cost  function  values 
between  the  minimax  design  and  the  optimal  fixed  design.  Referring  to  Table  5.1,  we  see 
that  the  optimal  cost  function  values  for  the  minimax  design  are  no  more  than  7%  larger 
than  those  for  the  optimal  fixed  design.  The  minimax  cost  function  value  and  damping 
coefficient  values  are  nearly  identical  with  those  obtained  for  fixed  optimal  design  for  the 
Puget  Sound  earthquake.  Note  that  for  the  El  Centro  earthquake,  although  the  minimax 
design  used  no  knowledge  of  the  El  Centro  ground  motions,  the  optimal  cost  function 
value  for  the  minimax  design  stmcture  is  within  2%  of  that  of  the  optimal  fixed  design, 
which  used  full  ground  motion  knowledge. 

•  peak  interstory  drifts.  Referring  to  Table  5.2,  we  see  that  for  the  San  Fernando  and 
Western  Washington  earthquakes,  the  minimax  design  exhibits  up  to  25%  larger  peak 
interstory  drifts  than  the  optimal  fixed  designs.  For  the  El  Centro  earthquake,  the 
minimax  design  exhibits  no  more  than  10%  larger  interstory  drifts.  For  the  Kem  County 
and  Puget  Sound  Earthquakes,  neither  of  the  designs  has  a  clear  advantage. 

•  peak  accelerations.  Referring  to  Table  5.3,  we  see  that  for  the  El  Centro  earthquake, 
the  minimax  design  sustained  no  more  than  10%  larger  peak  accelerations  than  the 
optimal  fixed  design.  For  the  Kem  County  and  Puget  Sound  quakes,  the  peak  accelera¬ 
tions  in  the  two  designs  are  nearly  equal,  while  for  the  San  Fernando  and  Western  Wash¬ 
ington  quakes,  the  optimal  fixed  designs  sustain  up  to  20%  more  peak  accelerations. 

•  average  accelerations.  Referring  to  Table  5.4,  we  see  that  for  the  El  Centro  and  Kem 
County  earthquakes,  the  minimax  design  sustained  no  more  than  8%  more  average 
accelerations  than  the  optimal  fixed  designs;  for  the  Puget  Sound  earthquake,  the  average 
accelerations  are  nearly  identical  for  the  two  designs;  and  for  the  San  Fernando  and 
Western  Washington  earthquakes,  the  optimal  fixed  designs  exhibit  up  to  15%  more 
average  accelerations  than  the  minimax  design. 

•  variable  damper  energy  dissipation.  Referring  to  Table  5.5,  we  see  that  for  the  El 
Centro,  Kem  County,  and  Puget  Sound  earthquakes,  neither  the  minimax  nor  fixed 
optimal  designs  have  a  clear  advantage,  while  for  the  San  Fernando  and  Western  Wash¬ 
ington  earthquakes,  the  minimax  design  requires  up  to  45%  more  energy  to  be  dissipated 
in  the  dampers  than  for  the  optimal  fixed  designs. 

•  main  structure  shear  force.  Referring  to  Table  5.6,  we  see  that  for  the  El  Centro  and 
Western  Washington  earthquakes,  the  minimax  design  sustained  up  to  18%  more  shear 
force  than  the  optimal  fixed  designs,  while  for  the  remaining  earthquakes,  neither  of  the 
designs  exhibits  a  clear  advantage. 


-31- 


5.  Evaluation  of  Numerical  Results 


•  peak  damper  force.  Finally,  referring  to  Table  5.7,  we  see  that  for  the  El  Centro  and 
Kern  County  earthquakes,  the  minimax  design  requires  no  more  than  11%  larger  damper 
forces  than  the  optimal  fixed  designs.  For  the  Puget  Sound  earthquake,  the  peak  damper 
forces  are  equal  for  both  designs;  and  for  the  San  Fernando  and  Western  Washington 
earthquakes,  the  optimal  fixed  designs  require  up  to  20%  larger  forces  to  be  applied  by 
the  variable  dampers  than  for  the  minimax  design. 

•  conclusions.  The  above  results  show  that  the  minimax  design  is  surprisingly  noncon¬ 
servative,  and  that  a  minimax  design  carried  out  using  a  small  number  of  earthquakes 
results  in  excellent  performance  on  earthquakes  not  used  in  the  design. 

evaluation  of  the  feedback  control  laws.  Next,  we  turn  to  an  evaluation  of  the 
performance  of  our  moving  horizon  feedback  laws  presented  in  Section  4,  so  as  to  deter¬ 
mine  how  well  they  perform  relative  to  one  another  and  to  the  performance  bounds  pro¬ 
vided  by  the  minimax  designed  stmcture  and  the  optimally  controlled  structure. 

•  cost  gap.  For  each  of  our  simulations  of  the  controlled  structure,  we  computed  cost 
function  values  using  the  same  cost  function  f(u)  (3.2.7)  that  was  used  in  the  optimal 
control  experiments,  but  with  u  equal  to  the  control  constructed  by  the  moving  horizon 
control  laws.  Referring  to  Table  5.1,  we  see  that  the  the  cost  function  values  for  the 
sampled-data  moving  horizon  law  simulations  range  from  1.3  to  2.7  times  as  large  as  the 
corresponding  cost  function  values  for  the  continuous  moving  horizon  law  simulations. 

With  the  exception  of  the  Puget  Sound  earthquake,  the  cost  function  values  for  the 
moving  horizon  controlled  structure  are  roughly  half  way  between  those  of  the  optimally 
controlled  stmcture  and  those  of  the  minimax  designed  stmcture.  For  the  Puget  Sound 
earthquake,  the  cost  function  value  for  the  moving  horizon  controlled  stmcture  is  6% 
larger  than  that  of  the  minimax  design. 

•  peak  interstory  drifts.  Referring  to  Table  5.8,  we  see  that  with  the  exception  of  x}, 
the  peak  drift  between  floors  2  and  3,  corresponding  to  the  El  Centro  earthquake,  the 
peak  interstory  drifts  of  the  stmcture  with  the  sampled-data  moving  horizon  control  law 
are  up  to  2.6  times  larger  than  those  of  the  stmcture  with  the  continuous  moving  horizon 
control  law. 

Comparing  the  interstory  drifts  of  the  stmcture  with  the  continuous  moving  horizon 
control  law,  in  Table  5.8,  with  the  interstory  drifts  of  the  minimax  design,  shown  in 
Table  5.2,  we  see  that  for  the  El  Centro  and  San  Fernando  earthquakes,  the  minimax 
design  exhibits  up  to  67%  larger  interstory  drifts  than  the  stmcture  with  the  continuous 
moving  horizon  control  law.  For  the  remaining  earthquakes,  the  stmcture  with  the  con¬ 
tinuous  moving  horizon  control  law  exhibits  no  more  than  51%  larger  interstory  drifts 
than  the  minimax  design.  The  use  of  optimal  control  almost  invariably  results  in  up  to 
40%  smaller  interstory  drifts  than  the  use  of  the  continuous  moving  horizon  control  law. 
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S-D  Moving  Horizon 

Moving  Horizon 

Quake 

Xd 

r  2 

Xd 

r3 

Xd 

Xd 

r2 

Xd 

r3 

Xd 

EC 

1.34 

1.46 

2.00 

0.78 

1.59 

1.31 

KC 

2.86 

3.17 

0.89 

1.50 

1.18 

PS 

1.89 

1.49 

0.77 

1.47 

1.20 

SF 

1.65 

1.72 

0.71 

1.27 

0.97 

WW 

1.93 

2.27 

1.19 

0.75 

1.58 

0.87 

Table  5.8.  Peak  Interstory  Drifts  (cm ) 

•  peak  accelerations.  Referring  to  Table  5.9,  we  see  that  except  for  the  ground  floors  of 
the  structures  in  the  Kern  County  and  Western  Washington  simulations,  the  peak 
accelerations  sustained  by  the  structure  controlled  by  the  moving  horizon  sampled-data 
control  law  are  up  to  15%  larger  than  the  corresponding  accelerations  sustained  by  the 
structure  controlled  by  the  continuous  moving  horizon  control  law. 

Referring  to  Tables  5.2  and  5.9,  we  see  that  for  the  El  Centro  earthquake,  the 
minimax  designed  structure  sustained  10%-42%  larger  peak  accelerations  than  the  con¬ 
tinuous  moving  horizon  controlled  structure.  For  the  Kern  County  and  Puget  Sound 
earthquakes,  neither  the  minimax  designed  structure  nor  the  continuous  moving  horizon 
controlled  structure  exhibits  a  clear  advantage.  For  the  San  Fernando  and  Western 
Washington  earthquakes,  the  continuous  moving  horizon  controlled  stmcture  sustained 
up  to  88%  higher  peak  accelerations  than  the  minimax  designed  structure. 

With  the  exception  of  the  Western  Washington  earthquake,  the  continuous  moving 
horizon  control  law  resulted  in  no  more  than  20%  higher  peak  accelerations  on  all  floors, 
than  the  optimal  control  law.  But  for  the  Western  Washington  earthquake  the  peak 
acceleration,  resulting  from  the  use  of  the  continuous  moving  horizon  control  law,  on  the 
ground  floor  was  higher  than  with  the  optimal  control  law,  while  the  accelerations  on  the 
upper  two  floors  were  lower  than  with  the  optimal  control  law. 


•  average  accelerations.  Referring  to  Table  5.10,  we  see  that  the  sampled-data  moving 
horizon  controlled  structure  sustained  up  to  70%  larger  average  accelerations  than  the 
continuous  moving  horizon  controlled  structure,  for  all  earthquakes  considered,  except 
Western  Washington,  where  neither  control  law  shows  a  clear  advantage. 

Comparing  the  structure  with  the  continuous  moving  horizon  control  law  and  the 
minimax  designed  structure,  see  Table  5.4,  we  see  that  for  all  earthquakes  other  than  the 
San  Fernando  earthquake,  the  continuous  moving  horizon  controlled  structure  sustained 
lower  average  accelerations  than  the  minimax  designed  structure  on  the  top  two  floors, 
but  higher  average  accelerations  on  the  ground  floor.  The  minimax  design  sustained  up 
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S-D  Moving  Horizon 

Moving  Horizon 

Quake 

-1 

^peak 

^peak 

r-3 

■^peak 

^peak 

^peak 

-3 

^peak 

EC 

579 

385 

392 

351 

279 

329 

KC 

460 

395 

111 

390 

317 

321 

PS 

452 

363 

364 

258 

221 

329 

SF 

417 

342 

453 

369 

341 

359 

WW 

283 

230 

229 

276 

226 

423 

Table  5.9.  Peak  Accelerations  {cmlsec^) 


to  40%  larger  accelerations  on  the  top  two  floors,  while  the  controlled  structure  sustained 
no  more  than  20%  larger  accelerations  on  the  ground  floor. 

Except  for  Western  Washington,  the  continuous  moving  horizon  controlled  struc¬ 
ture  sustained  up  to  26%  larger  accelerations  on  the  top  two  floors,  than  the  optimally 
controlled  stmcture,  while  the  optimally  controlled  structure  sustained  up  to  12%  higher 
accelerations  on  the  ground  floor.  For  the  Western  Washington  earthquake,  the  two  con¬ 
trol  laws  produced  almost  equal  average  accelerations. 


S-D  Moving  Horizon 

Moving  Horizon 

Quake 

-1 

^avg 

^avg 

••  1 
^avg 

^avg 

-3 

^avg 

EC 

102.9 

91.3 

94.3 

70.4 

64.5 

74.8 

KC 

70.2 

55.3 

46.0 

41.3 

36.5 

45.0 

PS 

61.6 

65.8 

80.2 

55.3 

53.2 

72.2 

SF 

59.3 

54.1 

80.8 

57.3 

50.4 

59.4 

WW 

38.1 

33.9 

40.5 

32.7 

34.5 

44.7 

Table  5.10.  Average  Accelerations  (cm /sec 

•  variable  damper  energy  dissipation.  Referring  to  Table  5.11,  we  see  that  the  continu¬ 
ous  moving  horizon  controlled  structure  dissipates  up  to  3;5  times  more  energy  in  each  of 
the  variable  dampers  than  the  sampled-data  moving  horizon  controlled  structure,  except 
in  the  Western  Washington  earthquake  simulations,  where  neither  control  law  has  a  clear 
advantage. 

Referring  to  Table  5.5,  we  see  that  the  minimax  designed  structure  dissipates  up  to 
2  times  more  energy  in  each  of  the  variable  dampers  than  the  continuous  moving  horizon 
controlled  structure,  except  for  the  second  floor  of  the  San  Fernando  simulations,  where 
the  continuous  moving  horizon  controlled  structure  dissipates  about  20%  more  energy 
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than  the  minimax  design. 

The  energy  dissipation  in  each  of  the  variable  dampers  of  the  moving  horizon  con¬ 
trolled  structtire  is  up  to  85%  larger  than  in  the  optimally  controlled  structure. 


S-D  Moving  Horizon 

Moving  Horizon 

Quake 

1^1 

E? 

Et 

Et^ 

EC 

0.83 

1.59 

2.25 

1.32 

4.06 

4.23 

KC 

0.65 

1.46 

1.62 

0.69 

2.10 

1.68 

PS 

0.25 

0.72 

1.15 

0.63 

2.00 

1.34 

SF 

0.27 

0.81 

0.84 

0.74 

2.80 

1.39 

WW 

0.50 

0.89 

0.76 

0.32 

1.04 

0.57 

Table  5.1 1.  Comparison  of  Total  Energy  Dissipation  in  Dampers  (iVm  x  10^) 

•  main  structure  shear  force.  Referring  to  Table  5. 12,  we  see  that  except  for  the  second 
floor  in  the  El  Centro  simulations,  the  use  of  the  sample-data  moving  horizon  control  law 
results  in  up  to  2.7  times  higher  shear  forces  than  the  use  of  the  continuous  moving  hor¬ 
izon  control  law. 

Comparing  the  peak  shear  forces  in  the  structure  with  the  continuous  moving  hor¬ 
izon  control  law  with  those  in  the  minimax  designed  structure,  shown  in  Table  5.6,  we 
see  that  for  the  El  Centro  and  San  Fernando  earthquakes,  the  minimax  designed  structure 
exhibits  up  to  73%  larger  shear  forces.  For  the  Kern  County  earthquake,  neither  structure 
has  a  clear  advantage;  and  for  the  Puget  Sound  and  Western  Washington  earthquakes,  the 
continuous  moving  horizon  controlled  structure  exhibits  no  more  than  55%  larger  shear 
forces  than  the  minimax  design. 

The  structure  with  the  continuous  moving  horizon  control  law  exhibited  up  75% 
larger  peak  shear  forces  than  the  optimally  controlled  stmcture  (see  Table  5.6). 


•  peak  damper  force.  Finally,  referring  to  Table  5.13,  we  see  that  for  all  but  the  San 
Fernando  earthquake,  the  sampled-data  moving  horizon  controlled  structure  exhibits  up 
to  66%  larger  forces  on  the  variable  dampers  than  the  continuous  moving  horizon  con¬ 
trolled  structore.  For  the  San  Fernando  earthquake,  neither  control  law  clearly  exhibits 
an  advantage. 

Comparing  the  continuous  moving  horizon  controlled  stmcture  with  the  minimax 
designed  stmcture,  see  Table  5.7,  we  find  that  for  the  El  Centro  earthquake,  the  dampers 
in  the  minimax  design  apply  up  to  50%  larger  forces  than  the  dampers  of  the  continuous 
moving  horizon  controlled  stmcture.  For  the  remaining  earthquakes,  the  dampers  of  the 
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S-D  Moving  Horizon 

Moving  Horizon 

Quake 

n 

B 

^91 

B 

EC 

IHBH 

2.9 

4.9 

1.5 

3.1 

3.2 

KC 

IBB 

5.6 

7.8 

1.8 

2.9 

2.9 

PS 

3.1 

3.7 

3.7 

1.5 

2.9 

3.0 

SF 

1.5 

3.2 

4.2 

1.4 

2.5 

2.4 

WW 

3.8 

4.5 

2.9 

1.5 

3.1 

2.1 

Table  5.12.  Comparison  of  Peak  Shear  Force  (N  x  10^) 


continuous  moving  horizon  controlled  stracture  apply  no  more  than  50%  larger  forces 
than  the  dampers  of  the  minimax  designed  structure. 

The  continuous  moving  horizon  controlled  structure  exhibited  no  more  than  42% 
higher  peak  forces  on  the  variable  dampers  than  the  optimally  controlled  structure. 


S-D  Moving  Horizon 

Moving  Horizon 

Quake 

B 

B 

B1 

B 

EC 

9.8 

13.7 

13.6 

5.9 

8.9 

9.3 

KC 

9.5 

14.9 

18.5 

lA 

11.5 

12.7 

PS 

7.4 

12.2 

13.3 

4.8 

7.8 

10.0 

SF 

7.3 

9.7 

12.0 

6.6 

10.7 

12.1 

WW 

7.4 

10.4 

10.0 

5.3 

7.9 

8.3 

Table  5.13.  Comparison  of  Peak  Damping  Force  (N  x  10^) 

•  conclusions.  Because  of  the  rather  small  gap  between  the  the  performance  gap 
between  the  realistic,  minimax  designed  stmcture  and  the  conceptual,  optimally  con¬ 
trolled  stracture,  finding  a  feedback  law  that  is  better  than  the  minimax  design  is  quite 
difficult.  Moving  horizon  laws  using  estimates  of  ground  motion  one  short  horizon 
ahead  seem  to  have  the  best  chance  of  falling  into  the  performance  gap.  However,  even 
using  unrealistic  projections  of  the  the  ground  motion,  our  simple,  sampled-data  moving 
horizon  feedback  law  led  to  worse  results  than  our  minimax  design.  On  the  other  hand, 
our  continuous  moving  horizon  feedback  law,  again  using  unrealistic  projections  of  the 
the  ground  motion,  resulted  in  better  performance  than  the  minimax  design,  particularly 
on  the  severe,  El  Centro  earthquake. 
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validation  of  “short  time”  experiments.  To  establish  confidence  in  our  experi¬ 
mental  results  that  examined  only  the  first  5  seconds  of  the  response  of  a  structure  to  a 
ground  motion  excitation,  we  performed  one  set  of  experiments  using  the  the  first  20 
seconds  of  the  El  Centro  earthquake.  The  experiments  performed  were  an  optimal  control 
computation,  a  simulation  of  the  response  of  the  minimax  designed  structure,  and  a  simu¬ 
lation  of  the  response  of  the  structure  controlled  using  the  continuous  moving  horizon 
control  law  (4.2).  Other  than  using  a  final  time  equal  to  20  seconds,  the  parameters 
used  for  the  optimal  control  experiment  were  the  same  as  those  used  in  Section  3.2,  and 
the  parameters  for  the  moving  horizon  simulation  were  the  same  as  those  used  in  Section 
4.  The  resulting  interstory  drifts,  floor  accelerations,  variable  damping  coefficients  are 
shown  in  Figures  5.6  -  5.13. 

Referring  to  Table  5.2,  Table  5.8,  and  Table  5.14,  we  see  that  the  peak  interstory 
drifts  in  the  20  second  experiments  are  very  close  to  those  obtained  in  the  5  second 
experiments. 


Xd 

Xd 

Xd 

Optimal  Control 

0.14 

1.33 

0.99 

Minimax  Design 

1.31 

1.91 

1.77 

Moving  Horizon 

0.83 

1.61 

1.17 

Table  5.14.  Comparison  of  Peak  Interstory  Drifts  (cm) 

Comparing  the  results  in  Table  5.15  with  with  the  corresponding  results  in  Table 
5.3,  we  see  that  the  peak  accelerations  for  the  20  second  and  5  second  optimal  control 
experiments  are  nearly  equal,  as  are  the  peak  accelerations  for  the  20  second  and  5 
second  simulations  using  the  minimax  designed  structure.  Referring  to  Table  5.9  and 
Table  5.15,  we  see  that  the  peak  accelerations  sustained  by  the  continuous  moving  hor¬ 
izon  controlled  structure  are  somewhat  larger  for  the  20  second  simulation  than  for  the  5 
second  simulation,  but  the  difference  is  not  large  enough  to  change  any  of  the  conclu¬ 
sions  obtained  from  the  5  second  simulation  results. 


f  ^ 

"^peak 

i'2 

^peak 

‘^peak 

Optimal  Control 

283 

253 

273 

Minimax  Design 

498 

417 

365 

Moving  Horizon 

375 

373 

337 

Table  5.15.  Comparison  of  Peak  Accelerations  (cmlsec^) 

Referring  to  Table  5.4,  Table  5.10,  and  Table  5.16,  we  see  that  the  average 
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accelerations  for  the  20  second  experiments  are  all  smaller  than  for  the  5  second  experi¬ 
ments,  but  the  conclusions  obtained  from  the  5  second  experiment  about  the  relative  size 
of  the  average  accelerations  also  hold  for  the  20  second  experiment  results. 


-  1 

X 

•^avg 

-2 

^avg 

-3 

^avg 

Optimal  Control 

40.4 

37.5 

47.2 

Minimax  Design 

61.0 

50.2 

46.3 

Moving  Horizon 

48.5 

42.9 

45.8 

Table  5.16.  Comparison  of  Average  Accelerations  {cmtsec^) 


We  conclude  from  the  above  results,  that  the  5  second  variable  structure  experi¬ 
ments  give  an  adequate  indication  of  the  performance  of  the  various  designs  over  a  20 
second  horizon.  Thus,  it  seems  quite  likely  that  the  conclusions  resulting  from  the  5 
second  experiments  will  be  equally  valid  for  longer  time  horizons. 

6.  Conclusion 

Using  the  responses  to  seismic  excitation  of  an  optimally  controlled  variable  struc¬ 
ture  and  of  a  minimax  optimally  designed  fixed  structure,  we  have  established  an  upper 
bound  on  the  achievable  performance  and  a  lower  bound  on  the  acceptability  of  a  control 
system  for  a  variable  damping  structure.  The  gap  between  the  upper  and  lower  bounds  is 
rather  small,  which  indicates  that  to  design  a  control  system  that  results  in  a  variable 
structure  that  is  clearly  superior  to  a  minimax  optimal  designed  fixed  structure  is  a  very 
difficult  task  indeed.  Our  experimental  results  indicate  that  a  controlled  variable  struc¬ 
ture  is  likely  to  perform  better  than  a  fixed  structure  in  the  case  of  moderate  to  severe 
earthquakes.  Further  reflection  leads  us  to  believe  that  controlled  variable  stractures  are 
likely  to  perform  best  at  sites,  such  as  landfills  and  dry  lake  beds,  where  resonances  can 
be  expected,  but  the  resonance  frequency  cannot  be  estimated  in  advance.  The  reason  for 
this  is  that  a  variable  structure  will  adjust  itself  for  an  impedance  mismatch  (i.e.,  for 
anti-resonance),  while  a  fixed  structure  is  not  capable  of  such  adaptation. 

In  addition,  we  have  found  that  even  fairly  simple  minimax  formulations  of  the 
design  of  a  fixed  structure,  based  on  a  small  number  of  ground  motions,  produce  very 
good  results,  in  the  sense  that  the  resulting  structure  meets  specifications  on  interstory 
drifts,  shear  forces,  and  floor  accelerations  not  only  for  the  earthquakes  considered  in  the 
design  process,  but  also  for  other  earfliquakes  scaled  to  comparable  intensity.  Sophisti¬ 
cated  minimax  design  techniques  (see,  e.g.,  [Pol.  1,2]),  using  direct  measures  of  perfor¬ 
mance  requirements  in  the  form  of  bounds  on  instantaneous  interstory  drifts,  shear 
forces,  accelerations,  etc.,  are  likely  to  produce  even  better  results. 
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Finally,  we  have  found  that  a  continuous  moving  horizon  control  law,  using  a  hor¬ 
izon  of  only  0.2  seconds,  produces  quite  acceptable  performance  in  a  variable  structure. 
However,  the  implementation  of  this  control  law  will  require  the  resolution  of  two  issues. 
The  first  is  that  of  speed  of  computation,  and  will  require  finding  suitable  hardware  as 
well  as  developing  particularly  efficient  software.  The  second  issue  is  that  of  ground 
motion  prediction  over  an  interval  of  0.2  seconds.  It  will  be  necessary  to  determine 
whether  such  a  prediction  can  be  made  by  using  ground  motion  monitoring  sensors 
located  a  small  distance  away  from  the  site. 
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Appendix:  Discretized  Gradient  Computation 

In  this  appendix,  we  will  derive  in  detail  the  formula  for  the  computation  of  the  gra¬ 
dient  of  the  cost  function  f  (u)  (3.2.4)  for  the  discretized  optimal  control  problem 
described  in  section  3.2.  We  begin  by  considering  the  more  general  problem  of  the  com¬ 
putation  of  the  gradient  of  the  cost  function  for  an  optimal  control  problem  with  dynam¬ 
ics  of  the  form 

x(r)  =  ^(x(r),n(r)),r6[0,ry.],x(0)  =  0,  (A.l) 

where  p.  6  M  with  M  defined  in  (3.2.1a),  and  x  e  R" .  The  cost  function  is  assumed  to 
be  of  the  form 

F(ii)  =  g(xHTf)),  (A.2) 

where  g  (•)  is  continuously  differentiable  and  x^(  )  is  the  solution  of  (A.l)  subject  to  the 
control  p.  We  assume  that  the  function  h  (*,  •)  is  continuously  differentiable  and  that 
•)  and  ft^(-,  •)  are  Lipschitz  continuous.  The  problem  to  be  considered  is 

p  minF(p).  (A.3) 

fieM  '  ' 

In  order  to  obtain  a  numerical  solution  to  P^,  this  problem  must  be  replaced  by  a 
consistent  finite  dimensional  approximation  [Pol.4].  To  obtain  such  an  approximation, 
we  discretize  the  time  interval  and  replace  the  control  by  the  finite  parametric  representa¬ 
tion  (3.2.2a,b)  as  shown  in  Section  3.2.  Once  again,  the  backwards  Euler  method  is  used 
to  integrate  the  system  equations: 

X  ((it-i-l)A)  =  X  (kA)  +  AK  (X  ((k+l)A),  u  (()t-i-l)A)),  x~  (0)  =  0,  (A.4) 

where  we  assume  that  the  equations  may  be  solved  exactly. 

The  discretized  cost  function  then  becomes 

F(u)  =  g(x-'^(NA)),  (A.5) 

where  {x  “(^A) )  represents  the  solution  of  the  difference  equation  (A.5)  subject  to 
the  control  sequence  u  =  {u(k)} 

To  compute  the  gradient  of  the  cost  function,  VF  (m  )  we  first  look  at  the  differential 
of  the  solution  to  the  difference  equation  (A.4)  at  time  N  A ,  Dx  “  (N  A  ;bu)  =  6x  (N  A), 
where  fork  =  0, 1 , . . .,  ,  5x  ^“  (^  A)  is  the  solution  of  the  following  difference  equation: 

6x~ (ik+l)A)  -  5x~  (kA)  =  AH^ix  “  ((A:-i-l)A),  u (;t+l))5x" i(k+l)A) 

-I-  AF  „  (x  “  ((^-i-l)A),  u  (^-i-1))5m  (^-1-1),  (A.6a) 

subject  to  5x  (0)  =  0.  Rearranging  terms,  we  get 

bx  ((^+1)A)  =  A;t+i5x  (kA)  +  (^+1),  (A.6b) 

where  =  (/  -  AFj^(x“((fc+l)A),M(/:-i-l))^^  andF;t+i  =  AF„(x  “((^+1)A),m(^-I-1)). 
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Let  C>“  (•,  •)  be  the  state  transition  matrix  for  (A.6b),  i.e.,  for  any  N  >  k  >  j  >0, 
0“(A:A,7  A)  is  the  solution  to  the  matrix  difference  equation 

<!>i(k+l)^jA)  =  At^lmA,jA),  <i>(jA,jA)  =  I.  (A.6c) 

From  (A.6b,c),  we  obtain  the  following  expression  for  dx  (A^  A); 

5f  (NA)  =  (^+1)A)A^+i5^+i6m  (k+l).  (A.6d) 

k=0 

Next,  applying  the  chain  rule  to  the  function  F (u)  =  g  (x  “ (N A)),  we  get 

DF  (M  ;  5m  )  =  ( (;c  “  (NA)),  Dx  “  (N A ;  5m  ) ).  (A.7a) 

Using  the  fact  that  Z)i'“(lVA;5M)  =  5jc ^“(NA),  and  noting  from  (A.6c),  that 
0(NA,  (jt:+l)A)A^^.i  =  ^(NA,kA),  we  see  that 

N-l  „ 

DF(m;5m)=  X  {B][+i(b(NA,kAf7^g(x'*(NA)),bu(k+l)) 

*=o 

=  Z  {B][+^'^(kA),8uik+l)),  (A.7b) 

fc=0 

where  for  k  =  0, 1, . .  .,N, 

p  “(ifcA)  =  ^{NA,kAfV^g(x‘‘(NA)).  (A.7c) 

The  computation  of  the  adjoint  p  “  (•)  by  computing  the  state  transition  matrix  0(-,  •) 

and  using  equation  (A.7c)  is  not  particularly  efficient,  and  in  fact  requires  n  times  more 
operations  than  are  required  by  the  more  efficient  scheme  described  below.  From  (A.7c) 
we  get 

p  “(A:A)-p  “((/:+l)A)  =  mNA.kAf  -<D(NA,()t+l)A)^]V^g(f  “(NA)).  (A.7d) 
Using  the  fact  that  <I>(N  A,  (A:+l)A)A;t+i  =  <b(N  A,  kA),  we  get  that 

0(N  A,  kA)-  0(N  A,  (it+l)A)  =  0(N  A,  kA)  AK  ^{x  “  ((it+l)A),  u  (A:+l)).  (A.7e) 
Combining  (A.7d)  and  (A.7e)  results  in 

p  “ {kA)-p  “ ((ik+UA)  =  AK^ix  “ i(k+l)A).  u (ik+l))^0(N A, k  A)^  (x  “ (NA)) 

=  AR^(x  “ ((^+1)A), M (/:+l))^ p  “ A),  (A.70 

subject  top  “(NA)  =  V^g(x“(NA)).  Thus,  we  see  that  the  adjoint  is  computed  back¬ 
wards  in  time,  again  using  a  backwards  Euler  formula. 

The  optimal  control  problem  (3.2.1)  considered  in  Section  3.2  has  an  integral,  rather 
than  endpoint  cost  function,  so  the  above  results  cannot  be  directly  applied  to  problem 
(3.2.1).  However,  as  we  will  show,  it  is  not  difficult  to  transform  an  optimal  control 
problem  with  an  integral  cost  into  a  problem  with  an  endpoint  cost.  Consider  the  system 
model  (2.1.2).  For  notational  convenience,  we  will  denote  the  system  equations  by 
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Hf)  =  h(z (t), z (0)  =  0,  (A.8a) 

and  the  cost  function  by 

e{z^it))dt,  (A.8b) 

where  z^{‘)  is  the  solution  of  the  differential  equation  (A.8a)  subject  to  the  control  4(  ). 
The  discretized  system  equations  are  then  given  by 

z ((ifc+l)A)  =  z(k A)  +  Ah  (zi{k+l)A),  u  ()t+l)),  z  (0)  =  0,  (A.9a) 

and  the  discretized  cost  function  is  given  by 

f{u)  =  A^e(J^{kA)),  (A.9b) 

Jfc=l 

where  { z“(fc  A) }  is  the  solution  of  the  difference  equation  (A.9a)  subject  to  the  con¬ 
trol  sequence  { ii  (^) } 

To  convert  this  to  a  problem  of  the  form  considered  in  (A.4,  A.5),  we  define  for 
k  =  0, 1, . .  the  augmented  state 

X  (kA)  =  (z(kA)J(k)),  (A.lOa) 

where  z(fcA)e  is  the  system  state  defined  by  equation  (A.9a),  and  /(^)G  R  is 
defined  by 

=  A2  e(z0'A)).  (A.10b) 

j=i 

The  augmented  system  equations  are  then 

F((A:+1)A) 
f(k+l) 

subject  to  X  (0)  =  0,  which  is  in  the  form  (A.4).  The  cost  function  /  (m  )  is  then  given  by 

f(u)={e^^,x'*iNA)),  (A.10d) 

where  e  R^^  denotes  the  13**  column  of  a  13x13  identity  matrix.  Thus,  we  see  that 
the  original  problem  has  been  converted  to  the  form  (A.4,  A.5). 

To  compute  the  gradient  for  the  discretized  optimal  control  problem  (A.9),  we  will 
first  apply  the  results  derived  for  the  problem  (A.4,  A.5)  to  the  augmented  system  (A.  10). 
From  (A.7b),  we  see  that 

Df  (u  ;du)  =  A'^  (^„(f  “((^+1)A),m(A:+1))V“(^A),5m(;:-i-1)).  (A.lla) 
k=o 

For  k  =  0, 1, . .  .,N,  let p‘*(kA)  equal  the  first  12  components  of  p^ikA),  and  let PgikA) 
be  the  13**  component  of p  “  (kA).  Then,  from  (A.70,  we  get 


z(lcA) 

+  A 

f(k) 

/i(F((ifc+l)A),M(it+l)) 

e(^((k+m)) 


(A.lOc) 
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p'^ikA) 

p:ikA)_ 

_p“((/:-H)A)_ 

+  A 


K  (r  ((^+l)A),  u  (k+l)f  V,e  (F“  (()k+l)A))^ 


0 


0 


p\kA) 

PeikA) 

(A.  11b) 


subject  to p  “  (A A)  =  (u)  =  e  From  the  above  equations,  we  can  conclude  that  for 

k  =  0, 1, . .  .,N ,Pg(kA.)  =  Pg(NA)  =  1.  Thus,  equation  (A. 11b)  becomes 


p‘*ikA)=Api(k+l)A)p''((k+l)A)  +  AAp((k+l)AW,e(z^((.k+l)A)f ,  (A.llc) 

subject  to p “ (N A)  =  0,  where  Ap  ((A:+1)A)  =  (/  -  Ah^  (z*  ((^+1)A),  u  (/:+l))^)“^ 

The  equation  for  Df  {u  ;  bu ),  may  be  similarly  simplified  by  substituting  the  expres¬ 
sion  for  ^  „  (f  “  ((^-t-l)A),'M  (^:+l))  into  equation  (A.l  la),  which  becomes 


Df  (u  ;bu)  =  aY  {\hu(^ ((/: +1  )A). u(k+l)f  0 
/:=0  -* 


p'^ikA) 

P^ikA) 


,bu(k+l)) 


=  (/i„(?“((^+1)A),m(/:+1))V(^A),5m(^+1)).  (A.lld) 

k=0 

Referring  to  the  finite  dimensional  scalar  product  formula  (3.2.6b)  and  equation  (A.lld), 
we  see  by  inspection  that  for  k  =  1,2,.  ..,N, 

Vf(u)(k)  =  h^(r(kAXu(k)f  p^iik-DA).  (A.12) 
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.1:  Floor  Accelerations  (alpha  =  0.0) 


Time  (seconds)  Time  (seconds) 


F.g.  2.3..t:  Floor  AccelCTMions  (alpha  =  0.«))  Fig.  3.1.1 :  Omund  Position  (EC) 


Fig.  3. 1 .3:  Ground  Acceleration  (EC)  _ Fig.  3.1. 5:  Scaled  Ground  Velocity  (KC) 


Time  (seconds) 


Fig.  3.1.7:  Sealed  Ground  Position  (PS)  Fig.  3.1.9:  Scaled  Ground  Acceleration  (PS) 
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Fig.  3.2.3:  Damping  Coefficients,  Optimal  Control  (EO  Fig.  3.2.5:  Shear  Force.  Optimal  Control  (EC) 


Time  (seconds) 


Fig.  3.3.1:  Inlerslofy  Drifts.  Minimax  Design  (EC)  ■»<<>*  Fig  3-3.3:  Energy  Dissipation  in  Dampers,  Minima  Design  (EC) 
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Fig.  3.3.5:  Damping  Force,  Minimax  Design  (EC)  ^'8-  Interstory  Drifts,  S-D  Moving  Horizon  (EC) 


Time  (seconds) 


^4,3:  Damping  Coefficient..  S-DMovingHorizon  TO  ,  Fig,  45:  Shear  Fcce.  S-D  Movinn  znz. 


Time  (seconds) 


Inleretoiy  Drifts,  Moving  Horizon  (EC)  *10’  Fig.  4,9:  Damping  Coefficients,  Moving  Horizon  (EQ 


nme  (seconds)  Time  (seconds) 


Fig.  4. 1 1 :  Shelf  Foice,  Moving  Horizon  (EC) _  _ Fig,  5. 1 :  Intenitofy  Drifts.  Optimal  Fixed  (EQ 


(N)»ioj 
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Fig.  5.3:  Energy  Dissipaiion  in  Dampers,  Optimal  Fixed  (EC)  Fig.  5.5:  Damping  Force.  Optimal  Fixed  (EC) 


Time  (seconds) 


Fig.  5.6:  Intcisiory  DriHs,  Optima)  Control  (EC)  *  10*^  Fig.  5.8:  Damping  Coefficients,  Optimal  Control  (EQ 
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